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1.0 Configuration Determination 


To determine the best configuration from all candidate configurations it was pessary 
first to calculate minimum system weights and performance To optimize the desgn 
neceaary to vary configuration-specific variables such as total system weight, thrust-to- 
weightTatiosJoum durations, fetal thrust available, and mass fraction for_ the' jystem. 
Ootimizine each of these variables at the same time is technically unfeasible and not 
necessarily mathematically possible. However, discrete sets rf datat can begenerai e 
which will eliminate many candidate configurations. From the P S 

remaining designs, a final configuration can be selected. 

To ease this optimization process an interactive spreadsheet wascreatedltacceptedinitial 
thrust to weight ratios, number of main engines, vehicle diy weight, and ^ al (p e ’ 
seDaration) bum time. From this information, it calculated maximum weights and 
propellant mass fractions available for each vehicle and the final post-separation u 

time for the orbiter. 

After selecting the number of engines and thrust to weight ratios for both the boosteir’“d *he 
orbiter the length of the initial bum was varied to maximize fuel available to the orbiter for 
the second phafebum. This maximum weight was then used to determine tank sizings 
which dictated total volume necessary within the booster and orbiter. 

Included are the three most important designs considered: one which closely approximates 
the design criteria set forth in a Marshall Space Flight Center study of the Shuttle II, the 
configuration used in the initial proposal; and the final configuration presented and 
analyzed in this report. Note that these three systems represent important iterative steps in 

our final design. 

A listing by cell of the formulas used to generate the aforementioned data is included for 
reference. 
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ARBITER DATA - 


Initial Boost Burn Time 

.Thrust/Weight 

1.2000 

Initial Boost Engines 

4.00 

Final Boost Burn Time 

Total Thrust at Takeoff 

1, 454,100.00 lbf 

Total Burn Time 

Available Takeoff Weight 

1,211,750.00 lbm 


Dry Weight 

180,000.00 lbm 

Tsl 

Payload 

60,000.00 lbm 

Tvac 

RCS Fuel (MMH/N04) 

20 , 000.00 lbm 

Wdot 

Total Orbiter Weight 

260, 000.00 lbm 

SG - Hydrogen 

Maximum Fuel Weight 

951,75 0.00 lbm 

SG - Oxygen 

Density of Water 

Fuel Available Initially 

951,750.00 lbm 

Max Fuel - Hydrogen 

Initial Boost Engines 

4.00 

Initial Fuel Consumed 

505,080.80 lbm 

Max Fuel - Oxygen 

Primary Fuel Available 

446, 669.20 lbm 

Volume - Hydrogen 

Auxilliary Fuel 

496, 95 2.57 lbm 

Volume - Oxygen 

Total Fuel Available 

94 3, 621.. 7 7 lbm 

Total Volume 

Final Boost Engines 

4.00 

Eq. Hydrogen Tanks 

Approximate Final Thrust 

1,740,000.00 lbf 

Eq. Oxygen Tanks 

Approximate Final T/W 

1.4456 

Minimum Length Necessary 

BOOSTER DATA - 



Thrust/Weight 

1.4000 


Initial Boost Engines 

5.00 


Total Thrust at Takeoff 

1, 817,625.00 lbf 


Available Takeoff Weight 

1,298,303.57 lbm 


Dry Weight 

160 , 000.00 lbm 

Max Fuel - Hydrogen 

RCS Fuel (MMH/N04) 

10,000.00 lbm 

Total Booster Weight 

170,000.00 lbm 

Max Fuel - Oxygen 

'Fuel Weight 

1,128,303.57 lbm 

Volume - Hydrogen 

V 


Volume - Oxygen 

V , * ' 

Fuel Available Initially 

1,128,303.57 lbm 

Total Volume 

Initial Boost Engines 

5.00 

Eq. Hydrogen Tanks 

Initial Boost Fuel 

631,351.00 lbm 

Eq. Oxygen Tanks 

Reservoir Fuel Pumped 

4 96, 95 2.57 lbm 

Minimum Length Necessary 

OVIRALX, DATA 


m 

Thrust/Weight 

1.3034 


Initial Boost Engines 

9.00 


Total Thrust at Takeoff 

3,271,725.00 lbf 


Available Takeoff Weight 

2,510,053.57 lbm 


Dry Weight 

340, 000.00 lbm 


Payload 

60,000.00 lbm 


RCS Fuel (MMH/N04) 

30,000.00 lbm 


Total Dry Weight 

430,000.00 lbm 


Fuel Weight 

2,080,053.57 lbm 


Fuel Available Initially 

2,080,053.57 lbm 


Initial Boost Engines 

9.00 


Initial Boost Fuel 

1,136,431.80 lbm 


Fuel Available Finally 

943,621.77 lbm 



133.00 sec 
248.48 sec 
3 81.48 sec 

363 , 525.00 lbf 
435 , 000.00 lbf 

949.40 lbm/sec 
0.0700 ft~3 
1.1490 ft A 3 
62.43 lbm/ft"3 

135,964.29 lbm 
815, 785.71 lbm 
31,112 . 40 ft*3 
11,372.68 ft'3 
42,485.07 ft-3 
88.0301 ft 
21.4203 ft 
109.4504 ft 


161,186.22 lbm 
967,117.35 lbm 
36,883.88 ft "3 
13,482.36 ft-3 
50,366.23 ft'3 
69.4704 ft 
25.3939 ft 
94.8643 ft 
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ORBITER DATA - 


Thrust/Weight 

1.2000 

Initial Boost Engines 

5.00 

Total Thrust at Takeoff 

1, 817,825.00 lbf 

Available Takeoff Weight 

1,514,687.50 lbm 

Dry Weight 

180 , 000.00 lbm 

Payload 

60 , 000.00 lbm 

RCS Fuel (MMH/N04) 

20,000.00 lbm 

Total Orbiter Weight 

260 , 000.00 lbm 

Maximum Fuel Weight 

1,254,687.50 lbm 

Fuel Available Initially 

1,254,687.50 lbm 

Initial Boost Engines 

5.00 

Initial Fuel Consumed 

564,883.00 lbm 

• 

Primary Fuel Available 

68 9, 794.50 lbm 

Auxilliary Fuel 

563,410.57 lbm 

Total Fuel Available 

1,253,205.07 lbm 

Final Boost Engines 

4.00 

Approximate Final Thrust 

1, 740,000.00 lbf 

Approximate Final T/W 

1.1499 

BOOSTER DATA - 


Thrust/Weight 

1.4000 

Initial Boost Engines 

5.00 

Total Thrust at Takeoff 

1,817,625.00 lbf 

Available Takeoff Weight 

1,298,303.57 lbm 

Dry Weight 

160,000.00 lbm 

RCS Fuel (MMH/N04) 

10,000.00 lbm 

Total Booster Weight 

170,000.00 lbm 

Fuel Weight 

1 , 128 ,303.57 lbm 

Fuel Available Initially 

1,128,303.57 lbm 

Initial Boost Engines 

5.00 

Initial Boost Fuel 

564,893.00 lbm 

Reservoir Fuel Pumped 

563,410.57 lbm 

OVERALL DATA - 


Thrust/Weight 

1.2923 

Initial Boost Engines 

10.00 

Total Thrust at Takeoff 

3,635,250.00 lbf 

Available Takeoff Weight 

2 , 812,991.07 lbm 

Dry Weight 

340,000.00 lbm 

Payload 

60,000.00 lbm 

RCS Fuel (MMH/N04) 

30,000.00 lbm 

Total Dry Weight 

430 , 000.00 lbm 

Fuel Weight 

2,382,991.07 lbm 

Fuel Available Initially 

2,382,991.07 lbm 

Initial Boost Engines 

10.00 

Initial Boost Fuel 

1,129,786.00 lbm 

Fuel Available Finally 

1,253,205.07 lbm 


Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


119.00 sec 
330.00 sec 
44 9.00 sec 

363 , 525.00 lbf 
435 , 000.00 lbf 

949.40 ibm/sec 
0.0700 ft ^3 
1.1490 ft*3 
62.43 lbm/ft~3 

179,241.07 lbm 
1 ,075,446.43 lbm 
41,015.32 ft A 3 
14 , 992.55 ft * 3 
56,007.87 ft*3 
116.0496 ft 
28.2383 ft 
144.2879 ft 


161, 186.22 lbm 
967,117.35 lbm 
36,883.88 ft~3 
13,482.36 ft~3 
50,366.23 ft^3 
69.4704 ft 
25.3939 ft 
94.8643 ft 
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ORBITER DATA - 


Thrust/Weight 

1.2000 

Initial Boost Engines 

4.00 

Total Thrust at Takeoff 

1, 454,100.00 lbf 

Available Takeoff Weight 

1,211,750.00 Ibm 

Dry Weight 

180,000.00 Ibm 

Payload 

60,000.00 Ibm 

RCS Fuel (MMH/N04) 

20,000.00 Ibm 

Total Orbiter Weight 

260, 000.00 Ibm 

Maximum Fuel Weight 

951,750.00 Ibm 

Fuel Available Initially 

951,750.00 Ibm 

Initial Boost Engines 

4.00 

Initial Fuel Consumed 

562,044.80 Ibm 

Primary Fuel Available 

389,705.20 Ibm 

Auxilliary Fuel 

544,897.09 Ibm 

Total Fuel Available 

934, 602.29 Ibm 

Final Boost Engines 

4.00 

Approximate Final Thrust 

1, 740,000.00 lbf 

Approximate Final T/W 

1.4544 

BOOSTER DATA ------- 


Thrust/Weight 

1.4000 

Initial Boost Engines 

6.00 

Total Thrust at Takeoff 

2,181,150.00 lbf 

Available Takeoff Weight 

1,557,944.29 Ibm 

Dry Weight 

140,000.00 Ibm 

RCS Fuel (MMH/N04) 

10,000.00 Ibm 

Total Booster Weight 

170,000.00 Ibm 

Fuel Weight 

1, 387,944.29 Ibm 

Fuel Available Initially 

1,387,944.29 Ibm 

Initial Boost Engines 

6.00 

Initial Boost Fuel 

843,047.20 Ibm 

Reservoir Fuel Pumped 

544,897.09 Ibm 

OVERALL DATA ------- 


Thrust/Weight 

1.3125 

Initial Boost Engines 

10.00 

Total Thrust at Takeoff 

3,435,250.00 lbf 

Available Takeoff Weight 

2,749,714.29 Ibm 

Dry Weight 

340,000.00 Ibm 

Payload 

40,000.00 Ibm 

RCS Fuel (MMH/N04) 

30,000.00 Ibm 

Total Dry Weight 

430,000.00 Ibm 

Fuel Weight 

2,339,714.29 Ibm 

Fuel Available Initially 

2,339,714.29 Ibm 

Initial Boost Engines 

10.00 

Initial Boost Fuel 

1, 405,112.00 Ibm 

Fuel Available Finally 

934,402.29 Ibm 


Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


148.00 sec 

246.10 sec 

394.10 sec 

363 , 525.00 lbf 
435 , 000.00 lbf 

949.40 lbm/sec 
0.0700 ft *3 
1.1490 ft *3 
62.43 lbm/ft-3 

135,964.29 Ibm 
815,785.71 ibm 
31,112.40 ft^3 
11,372.68 ft*3 
42,485.07 ft~3 
88.0301 ft 
21.4203 ft 
109.4504 ft 


198,280.61 Ibm 
1, 189, 683.67 Ibm 

45.372.10 ft*3 

16.585.10 ft~3 
61, 957.20 ft^3 

85.4579 ft 
31.237 9 ft 
116.6958 ft 
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ORBITER DATA - 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
Payload 

RCS Fuel (MMH/N04) 

Total Orbiter Weight 
Maximum Fuel Weight 


1.2 

5 

-B3*E6 

-B4/B2 

180000 

60000 

20000 

-B6+B7+B8 

-B5-B9 


Ibf 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 


Fuel Available Initially 
Initial Boost Engines 
Initial Fuel Consumed 
Primary Fuel Available 
Auxilliary Fuel 
Total Fuel Available 
Final Boost Engines 
Approximate Final Thrust 
Approximate Final T/W 

BOOSTER DATA - 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
RCS Fuel (MMH/N04) 

Total Booster Weight 
Fuel Weight 

Fuel Available Initially 
Initial Boost Engines 
Initial Boost Fuel 
Reservoir Fuel Pumped 

OVERALL DATA - 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
Payload 

RCS Fuel (MMH/N04) 

Total Dry Weight 
Fuel Weight 


-B10 

-B3 

-B13 *18*X2 
-B12-B14 

— B3 4 

—B15+B16 
4 

-B18*X7 

— B19/(B9+B17) 


1.4 

5 

—B2 3 * E 6 
-B24/B22 

160000 

10000 

-B26+B27 

-B25-B28 

— B2 9 

-B23 

— B3 2 *E8 *X2 
-B31-B33 

-(B4+B24)/(B5+B25) 

-B3+B23 

-B4+B24 

-B38/B36 

-B6+B26 

— B7 

-B8+B27 
—B40+B41+B42 
-B10+B29 


lbm 

lbm 

lbm 

lbm 

lbm 

lbf 


lbf 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 


lbf 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 


Fuel Available Initially 
Initial Boost Engines 
Initial Boost Fuel 
Fuel Available Finally 


-B12 +B31 
—B37 

— B4 7 *E2*E8 
-B46-B48 


lbm 

lbm 

lbm 
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Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


119 

sec 

-B17/(318*18) 

sec 

—12 >13 

sec 

363525 

ibf 

435000 

ibf 

949.4 

lbm/sec 

0.07 

ft A 3 

1.149 

ft *3 

62.43 

lbm/ft*3 

-B12/7 

lbm 

-B12 * 6/7 

lbm 

-113/(19*111) 

ft *3 

-114/ (110*111) 

f t A 3 

-115>116 

ft~3 

-115/ (PI ()/4*15 A 2)/2 

ft 

-116/(PI ()/4*2 6 A 2) 

ft 

—118>119 

ft 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


-B29/7 
«B29*6/7 
-127/(19*111) 

-£28/(£10*111) 
-I29+I30 

-129/(PI()/4*26 A 2) 
-130/(PI()/4 *2 6 A 2) 
—132 >133 


Ibm 
lbm 
ft "*3 
ft *3 
ft~3 
ft 
ft 
ft 
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2.0 STME Performance Calculation Summary 

A model of the total thrust available from the Space Transportation Main Engine 3 used on 
the Shuttle II was necessary to evaluate its launch profile. Thrust is a function of y 
SkrinclST^kexpansion ratio, exit velocity, ambient pressure and exit 

Stal ?** STME uses a translatable skirt, the expansion rata, of the engine could 

be set at discrete values of35:1,70:1, and 150:1. To insure that the^engine* ope 
efficiently during all points in the launch, it was necessaiyto calculate the altitudes at 
which tae skirt should be translated. This analysis was performed by calculating thrust 
coefficients for all three nozzles and finding the points of intersection. 

These calculations were based on the data presented in the Aeroject contractor's report for 
tae ^ ?962 standard atmosphere was used as a model ^pressured density as 

functions of altitude. Most of the engine data was solved using simple, iterative 
techniques, isentropic relations, and graphical methods. 
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2.1 Analysis 

Given an expansion ratio for a particular nozzle and an appropriate value for the 
coefficient of specific heats, the exit Mach number may be determined. 

Assuming y = 1.24 for a mixture of liquid oxygen and liquid hydrogen in a 6:1 mass ratio 
and the first (fixed) nozzle expansion ratio of 35:1, this isentropic relation holds: 

Aft _ _L- ( I f Ale 1 ”! z( |-i \ 

A* yJl J 

This function may be solved iteratively to obtain exit Mach numbers for the nozzles. This 
was solved using the included iterative, and the following were obtained as the exit Mach 

numbers: 


Nnmie Exp ansion Ratio 
35:1 
70:1 
150:1 


rerit Mach Number 
4.369 
4.931 
5.581 


Next, the flow in a nozzle may be assumed to be very nearly isentropic, such that the 
chamber pressure may be related to the nozzle exit pressure through the isentropic relation 

Given a chamb er pressure of 2385 psia under normal power levels and 3125 psia under 
emergency conditions, the following exit pressures were obtained. 


Pe (NFL) PefEPL) 

For e = 35:1 5.069 6.642 

For e = 70:1 2.058 2.696 

For e = 150:1 0.771 1.010 

For a given expansion ratio, a nozzle will have a maximum efficiency at a certain design 
altitude. At this altitude, the exit pressure is equal to the ambient pressure of the 
atmosphere. From the 1962 standard atmosphere, the design altitudes were determined to be 

hd (Normal Power Level) 

For e = 35:1 26,700 ft 

Fore = 70:1 45,800 ft 

Fore =150:1 66,300 ft 



20,400 ft 
40,200 ft 
60,600 ft 


In order to obtain the proper altitudes at which the sections of the nozzle should be extended, 
the following procedure is followed. 


The thrust coefficient is defined as: 



Since thrust varies with altutde, this expression may be rewritten as. 




£z} 
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For example, consider the nozzle with an expansion ratio of 70:1 at an altitude of 50,000 feet. 


c T * lL'*-* x **K* m +V ~ J 

Similarly, thrust coefficients for other nozzle and atmospheric conditions may ’be 
calculated. As shown in Figure 2.1.1, the points where these graphs orwiindicate 
altitudes where the nozzles should be translated for maximum efficie cy. 

For Normal Power Levels, the translation altitudes were found to ^36,900 feet and 56,500 
feet For Emergency Power Levels, the altitudes were 31,100 feet and 50,800 feet. 


Z.O60 ' i-OSZ. 


(-IO) 
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55 


29000 


30000 


31000 


32000 


33000 


34000 


35000 


36000 


37000 


38000 


39000 


40000 


41000 


42000 


43000 


44000 


45000 


46000 


47000 


48000 


49000 


50000 


51000 


52000 


53000 


54000 


55000 


B 


P 


14.7000 


14.1700 


13.6600 


13.1700 


12.6900 


12.2300 



NPL, £-35 


1.6185 


1.6263 


1.6338 


1.6410 


1.6480 


1.6548 


T"T 

6000 

11.7800 

1.6614 


ni 

11.3400 

1.6679 

in 

8000 

10.9200 

1.6740 


9000 

10.5000 

1.6802 

term 

10000 

10.1100 

1.6859 

HI 

11000 

9.7200 

1.6916 

mvm 

12000 

9.3460 

1.6971 

mrm 

13000 

8.9840 

1.7024 

in 

14000 

8.6330 

1.7076 

EH 

15000 

8.2940 

1.7126 


16000 

7.9650 

1.7174 

mrm 

17000 

7.6470 

1.7220 

if 

18000 

7.3390 

1.7266 

m 

19000 

7.04 l~ol 

1.7309 


20000 

6.7540 

1.73521 

23 

21000 

6.4750 

1.7392 

m 

22000 

6.2070 

1.7432 


23000 

5.9470 

1.7470 


24000 

5.6960 

1.7507 

WSM 

25000 

5.4540 

1.7542 


26000 

5.2200 

1.7577 

i 2 9 

27000 

4.9940 

1.7610 

eh 

28000 

4.7770 

1.7642 


4.5670 


4.3640 


4.1690 


3.9810 


3.8000 


3.6260 


3.4580 


3.2970 


3.1420 


2.9940 


2.8540 


2.7200 


2.5920 


2.4710 


2.3550 



2.1390 


2.0390 


1.9430 


1.8520 


1.7650 


1.6820 


1.6030 


1.5280 


1.4560 


1.3880 


1.3230 


1.7702 


1.7731 


1.7758 


1.7785 


1.7811 


1.7835 


1.7859 


1.7882 


1.7903 


1.7924 


1.7943 


1.7962 


1.7980 


1.7997 


1.8013 


1.8029 


1.8043 


1.8058 


1.8071 


1.8084 


1.8096 


1.8107 


1.8118 


1.8129 


1.8139 


1.8149 


EPL,E-35 


1.6696 


1.6756 


1.6813 


1.6868 


1.6921 


1.6973 


1.7023 


1.7073 


1.7120 


1.7167 


1.7210 


1.7254 


1.7296 


1.7336 


1.7376 


1.7414 


1.7451 


1.7486 


1.7521 


1.7554 


1.7586 


1.7617 


1.7647 


1.7677 


1.7705 


1.7732 


1.7758 


1.7783 


1.7808 


1.7831 


1.7854 


1.7876 


1.7897 


1.7917 


1.7937 


1.7955 


1.7973 


1.7991 


1.8007 


1.8023 


1.8038 


1.8052 


1.8066 


1.8079 


1.8091 


1.8103 


1.8114 


1.8125 


1.8135 


1.8145 


1.8154 


1.8163 


1.8172 


1.8180 


1.8187 


1.8194 


HP L, Z-7 


1.4493 


1.4649 


1.4798 


1.4942 


1.5083 


1.5218 


1.5350 


1.5479 


1.5602 


26 


1.5840 


1.5955 


1.6064 


1.6171 


1.6274 


1.6373 


1.6470 


1.6563 


1.6653 


1.6741 


1.6825 


1.6907 


1.6986 


1.7062 


1.7136 


1.7207 


1.7275 


1.7342 


1.7405 


1.7467 


1.7527 


1.7584 


1.7639 


1.7692 


1.7743 


1.7793 


1.7840 


1.7885 


1.7929 


1.7970 


1.8009 


1.8047 


1.8082 


1.8116 


1.8149 


1.8180 


1.8209 


1.8237 


1.8264 


1.8289 


1.8314 


1.8337 


1.8359 


1.8380 


1.8400 


1.8419 


ZPL,Z-70l HP 


1.5515 


1.5633 


1.5748 


1.5857 


1.5965 


1.6068 


1.6169 


1.6267 


1.6361 


1.6455 


1.6543 


1.6630 


1.6714 


1.6795 


1.6874 


1.6950 


1.7023 


1.7094 


1.7163 


1.7230 


1.7295 


1.7357 


1.7417 


1.7475 


1.7532 


1.7586 


1.7638 


1.7689 


1.7737 


1.7784 


1.7830 


1.7874 


1.7916 


1.7956 


1.7995 


1.8033 


1.8069 


1.8104 


1.8137 


1.8168 


1.8198 


1.8227 


1.8254 


1.8280 


1.8305 


1.8328 


1.8351 


1.8372 


1.8393 


1.8412 


1.8431 


1.8448 


1.8465 


1.8481 


1.8497 


1.8511 


L, E-1501 EPL 
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Figure 2.1.1: Plot of Thust Coefficients vs. Altitude (NPL) 
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Figure 2.1.2: Thust Coefficients vs. Altitude (EPL) 
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3.0 Mission Analysis 

To prove the flight worthiness and technical feasibility of the ASTS as proposed, a launch 
profile was constructed and simulated. This simulation consisited of a dynamic model of 
the vehicle’s instantaneous position, velocity, acceleration, thrust, weight, orientation, 
and flight path with respect to the horizon. In addition to this vehicle-specific data, 
atmospheric conditions such as temperature, pressure, and density were modeled. 


The basic format of this simulation was a fourth-order Runge-Kutta integration scheme. A 
number of asumptions were necessary to develop this model of the launch trajectory. First, 
the model is constrained by a run time of approximate seven minutes and a total distance 
downrange not exceeding 800 nautical miles, so the surface of the earth was modeled as 
being flat. In addition, rotational effects of the earth were neglected. This assumption 
cannot degrade the validity of results, since it could only create conditions more favorable 
to a due east launch from Kennedy Space Center. Drag forces were neglected for two 
principal reasons: they can be considered small in comparison to thrust forces necessary 
for a ballistic trajectory, and they would not exceed the advantages gained by 
incorporating the earth's rotation. 

All initial conditions were set based on values calculated from the design process outlined 
in Section 1.0. Bum times, initial weights, total thrust available, and fuel available are the 
the most noteworthy examples. Thrust was computed as a funciton of altitude, based on 
translation altitudes for the extendable skirt on the STME nozzles, pressure forces, and 
known values of sea level and vacuum values. It was estimated that the maximum 
acceptable pitch rate was 0.5 degrees per second to retain stability during launch. This is 
the value used in the program for all pitch manuevers. 

Before finalizing the launch profile, it was necessary to establish the range of orbits into 
which the orbiter could be successfully launched. Based on current STS mission analyses 
and estimations of maximum altitude and velocity available at this altitude, our launch 
goal was projected to be an altitude of 100 nautical miles and a velocity equal to that of the 
local circular speed. Inclination was varied during flight both to achieve this orbit and to 
place the booster in an alignment practical for separation. 









program AE447Prog; 

($MC68020+} 

{SMC68881+) 

{SSETC Elems881:-truel 


Written by Del Johnson (in TML Pascal II vl-00 under MPW v2. 02) _ __ __ 

Written””" AE447-AE449 Shuttle II design g roup (Fall 1988 - s P rir ^_j;f^_ 

;"rp''g"«'irdeII^d‘"“";r"r""hrof d,,ign configuration^for the 

Shuttle II Launch System. Design data was used in the model. ^® s 

execution of these programs was used in an iterative process to fine-tune the 

model as presented in the report. __ 

Definition of variables u[i]: 


{ Integral variables— } 


u ( 
u( 
u ( 
u [ 
u( 
u ( 


1 ] 

2 ] 

3] 

4] 

5] 

6 ] 


u[x] 
u (xd] 
u[y] 
u[yd] 
u [th] 
u [thd] 


Location along the x-axis 

Time rate of change of location along the 
Location along the y-axis 

Time rate of change of location along the 

Angle with respect to the y-axis 

Time rate of change of angle with respect 


x-axis 

y-axis 

to the y-axis 


( Other variables— } 
u[ ] - u[alpha]- Angle of attack 
u[ ] - u[gamma]- Flight path angle 

u[ ] - u[wt] - Total vehicle weight 

u[ ] « u[wtd] - Time rate of change of total vehicle weight 

u( ] - u[dens] - Density of the atmosphere 

u[ ] - u[press] - Density of the atmosphere 
u[ ] - u[ thrust] - Thrust produced by all engines firing 


*> 


uses MemTypes,SANE; 
{ 


const 


{ Number of equations- ) 

numlntEquations - 6; 

numTotEquations — 11; 


{ Time Considerations- } 


tbegin 

tend 

dt 

cint 

cols 


0.00 

420.00 

0.10 

10 

4 


( seconds } 
( seconds } 
{ seconds } 



Subscripts 


of the variables of motion— } 


I 
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xd ■ 1; 

x * 2; 

yd * 3; 

y *4; 

wt ■ 5; 

th - 6; 

vtot ■ 7; 

gamma * 8; 

thrust * 9; 

press * 10; 

dens * 11; 

Separation * 1 
Shutdown " 2 


StartPitchl ■ 3 
StopPitchl " 4 
StartPitch2 - 5 
StopPitch2 ■ 6 


{ Constants related to the problem- } 

{ Mathematical constants.- } 

conv - 57.29577951; { degrees/radians } 


{ 


Physical 

Re 

gO 

rho 

rhoO 

Psl 


constants.- } 

- 20938912.00000; { 

- 32.17400; { 

- 0.0023768; 

- 0.0023768; 

- 14.70000; { 


feet } 

feet/second A 2 } 

{ pound-mass/feet A 3 } 

{ pound-mass/feet''3 } 

pound-force/inches A 2 } 


{ Vehicle characteristics.- } 

thrO - 363535.00000; { Ibf } 

wtdO - -949.40000; { Ibm/sec } 


pressAlt * 25000 
dens Alt * 25000 

Ael - 3247 
Ae2 - 6503 
Ae3 - 14526 


00000; { feet } 
00000; { feet } 

22085; { inches A 2 } 
88219; { inches A 2 } 
72443; j inches A 2 } 


designAltlNPL 
de s ignAlt 2NP L 
de s ignAl 13NP L 
designAlt1EPL 
designAlt2EPL 
de s ignAlt 3 EPI* 


26700.00000; { feet } 
45800.00000; { feet } 
66300.00000; { feet } 
20400.00000; { feet } 
40200.00000; { feet } 
60600.00000; { feet } 


transAltlNPL - 36900.0000; { feet } 
transAlt2NPL - 56500.0000; { feet } 
transAltlEPI* - 31100.0000; { feet } 
transAlt2EPL - 50800.0000; { feet J 


PelNPL - 5.069; { lbf/feet / '2 } 
Pe2NPL - 2.058; { lbf/feet A 2 } 


4 
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Pe3NPL - 0.771 
PelEPL - 6.642 
Pe2EPL - 6.296 
Pe3EPL - 1.010 

thrMaxlNPL - 4 
thrMax2NPL - 4 
thrMax3NPL ” 4 

thrMaxlEPL * 5 
thrMax2EPL ■ 5 
thrMax3EPL “ 5 


1; { lbf/feet A 2 } 

2; { lbf/feet A 2 } 

6; { lbf/feet~2 } 
0; { lbf/feet A 2 } 

411259.15000; 

424577.04000; 

435000.00000; 

545789.15000; 

556312.08000; 

580000.00000; 


{ Problem intial conditions™ } 


xdO 

J 

0.00; 

( 

feet/second 

* 

xO 

- 

0.00; 

{ 

feet } 


ydO 

- 

0.00; 

{ 

feet/second 


yO 

- 

0.00; 

{ 

feet } 


thO 

- 

90.00; 

{ 

degrees } 


vtotO 

« 

0.00; 

{ 

feet/second 

} 


type 

uArray 
nameArray 
printrange 


array [1..numTotEquations] of extended; 
array [1..numTotEquations] of str255; 

1..32; 


var 


outflie 
tab 
3p 


char; 


array [1 • .nuniTotEquations] of integer; 


DryWtBooster, DryWtOrbiter,wtO : extended; 


extended; 


time 

k 


array [1..6] of extended; 

array [1.•4 r 1.•numlntEquations] of extended 


u,uu 
if j/m,n 
printset 


uArray; 

integer; 

set of printrange; 


{ Variables of motion^ } 

wtd, thr, thd : extended; 

thrMax : extended; 


function g(alt : extended): extended; 
begin 

g g0*sqr (Re)/sqr (Re+alt) ; 

end; 

function density(alt : extended): extended; 
begin 

density :* rho0*exp(—aIt/densAlt); 

end; 

function pressure(alt : extended): extended; 
begin 

pressure :- Psl*exp(-alt/pressAlt); 

end; 
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function pitchRate(t : extended): extended; 
begin 

if (t < Time[StartPitchl]) then pitchRate 0.00 else 

if (t < Time[StopPitchl] ) then pitchRate -0.50 else 

if (t < Time[StartPitch2]) then pitchRate :* 0.00 else 

if (t < Time[StopPitch2] ) then pitchRate -0.50 else 

if (t < tend) then pitchRate 0.00; 

end; 

function thrustAvail(alt : extended): extended; 
var 

Ae,thrl : extended; 

begin 

if (alt<transAltlNPL) then begin 
Ae Ael; 
thrl :» thrMaxlNPL; 

end 


if (alt<transAlt2NPL) then begin 
Ae Ae2; 
thrl :* thrMax2NPL; 

end 

else begin 

Ae :• Ae3; 

thrl :* thrMax3NPL; 

end; 

thrustAvail :** thrl - pressure(alt)*Ae; 

end; 

procedure initVariables; 
begin 

t :- tbegin; 
u[yd] ydO; 
u [y] yO; 

u[xd] xdO; 
u[x] xO; 

u[wt] :* wtO; 
u(th] :* thO; 
u[gamma] :* 90; 

u[vtot] sqrt(sqr(u[yd])+sqr(u[xd])); 
u[press] pressure(u[y]); 
u[dens] :» density(u[y]); 

end; 

I...................-.............................- 

procedure print(style : integer); 
var 

i : integer; 
begin 

case style of 
1 : begin 

write(outfile, f t ■ f ,t:6:2); 

for i :■ 1 to numTotEquations do begin 

if (i in printset) then write(outfile, 1 

-',u[i]:12:2); 

if (i in printset) then write( 1 

-*, u ti] : 12 :2) ; 

end; 


u[ 9 ,i:2, f ] 
u[ 9 ,i : 2 , 9 ] 


1 


} 


} 


1 
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'); 


writeln(outfile); 
writeln; 
end; 

2 : begin 

i :■ 0; 

writeln(outfile,'Time ■ ',t:7:2); 
writeln ( 'Time *■ ' ,t:7;2); 

for j ;■ 1 to numTotEquations do 
if (j in printset) then begin 
i : m i+1; 

write(outfile,'u[\ j:2,'] ,u( j]:12:2,' 

write( 'u[',j:2,«] -',u[j]:12:2,' 

if (i mod cols - 0) then begin 
i :« 0; 

writeln(outfile); 
writeln; 
end; 
end; 

if (i<>0) then begin 
writeln(outfile); 
writeln; 
end; 
end; 

3 : begin 

write(outfile,t:7:2,tab); 

for j :■ 1 to numTotEquations do 

if (j in printSet) then write(outfile,u(j]:12:sp(j],tab); 

writeln(outfile); 
end; 
otherwise; 
end; (case) 
end; 

{ Total function defining variables in terms of one another... ) 
function f(j: integer; t: extended; u: uArray): extended; 

begin 


j of 
xd 
x 

yd 

y 

wt 

th 


f 

f 

f 

f 

f 

f 


- thr*cos(u[th]/conv)*g(u[y])/u[wt]; 

» u[xd]; 

- thr*sin(u[th]/conv)*g(u[y])/u(wt] - g(u[y]); 

- u [yd] ; 

- wtd; 

- thd; 


otherwise; 


end; 


end; 


begin 


Textbook(nil) ; 
tab chr(09); 

for j :■ 1 to numTotEquations do 

sp[j3 :* 2; 

sp[press] :* 8; 
sp[dens] :* 8; 


open (outfile, 1 l/AE448prog.dat ); 
rewrite(outfile); 
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m numlntEquations; t Number of integration equations } 

n trunc((tend-tbegin)/dt); { Number of steps } 

{ Set names of all variables.- } 

printset [xd,x,yd, y,wt, vtot,th,gamma,thrust,press,dens]; 

{ Set design variables.- } 

wtO 2769714.29; { lbf } 

DryWtBooster :~ 170000.00; { lbf } 

DryWtOrbiter : m 260000.00; { lbf } 

Time [Separation] :*■ 148.00; { seconds } 

Time[Shutdown] :* 394.00; { seconds } 

Time[StartPitchl] 20.00; { seconds } 

Time[StopPitchl] 140.00; ( seconds } 

Time[StartPitch2] 160.00; { seconds ) 

Time[StopPitch2] :■ Time[StartPitch2] +38.00; [ seconds } 

( Set initial conditions of all variables- } 
initVariables; 

{ Echo initial conditions of all variables in the printing set— < 
print(3); 

{ Beginning of time-step loop— } 
for i 1 to n do begin 

{ Computation of other f non time—dependent variables— } 

if (t — Time[Separation]) then u[wt] : m u[wt] - DryWtBooster; 

thrMax :« thrustAvail(u[y]); 

if (t < Time [Separation]) then thr 10*thrMax else 

if (t < Time[Shutdown] ) then thr 4*thrMax else 

if (t < tend) then thr 0*thrMax; 

u[press] pressure(u[y]); 

u[dens] density[u[y]); 

u[thrust] ;■ thr; 

if (t < Time[Separation]) then wtd : m 10*wtd0 else 

if (t < Time[Shutdown] ) then wtd 4*wtd0 else 

if (t < tend) then wtd 0; 

thd pitchRate(t); 

if (u[th] > 360.00) then u[th] u[th] - 360.00; 

if (u[thj < -360.00) then u[th] :■ u[th] + 360.00; 

{ Runge-Kutta 4th-Order integration loop of all variables- ) 

for j ;■ 1 to a do 

k[l,j] dt*f(j,t,u); 

for j 1 to m do 

uu[ j] ;“ u[ j] + k[l,j]/2; 

for j :■ 1 to m do 

k[2,j] dt*f(j,t+dt/2,uu); 
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u[gamma] 


for j :» 1 to m do 

uu[j] :■ u[j] + k[2,j]/2; 

for j :« 1 to m do 

]c[3, j] dt*f (j,t+dt/2,uu) ; 

for j :« 1 to m do 

uu[j] :“ u[j] + k [ 3, j ]; 
for i 1 to m do 


k(4,j] :■ dt*f(j,t+dt,uu); 
for j :■ 1 to b do 

u [j] u[j] + (k [ 1 , j]+ 2 *k[ 2 ,j]+2*k [3/}]+k[4, 




t :■ tbegin + i*dt; 

u[vtot] :■ sqrt (sqr (u[xd])+sqr(u[yd])); 

if (u[xd] > 0) then u[gamma] :■ arctan (u[yd]/u[xd])*conv 

90.00; 


{ Printing parameters.- } 

if (i mod cint - 0) and (printset <> []) then print(3); 


end; 

close(outfile); 


end. 


else 
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l/AE448data 


Tima 

0.00 
1.00 
2.00 
3.00 
4,00 
5.00 
6.00 
7.00 
8.00 
9.00 
10.00 
11.00 
12.00 
13.00 
14,00 
15.00 
16.00 
17.00 


18.00 
19.00 
20.00 
21.00 
22.00 
23.00 
24.00 
25.00 
26.00 
27.00 
28.00 
29.00 
30.00 
31.00 
32.00 
33.00 
34.00 
35.00 
36.00 
37.00 
38.00 
39.00 
40.00 
41.00 
42.00 
43.00 


B 

C 

D 1 

E 


X-Vel 

X-Poa 

Y-Val 

Y-Pos 

Maiahtj 

0,00 

0.00 

0.00 

0.0 of 

2769714.2^] 

0.00 

0.00 

10.13 

5.05 

2760220.29) 

0.00 

0.00 

20.40 

20.30 

2750726.29? 

0.00 

0.00 

30.83 

45 • 91| 

2741232.29] 

0.00 

0.00 

41.41 

82. Olj 

2731738.29] 

0.00 

0.00 

52.15 

128.78[ 

2722244.29] 

0.00 

0.00 

63.05 

186.36) 

2712750.29] 

0.00 

0.00 

74.11 

254.93[ 

2703256.29] 

0.00 

0.00 

85.35 

334.65 

2693762.29] 

0.00 

0.00 

96.75 

425.68) 

2684268.29] 

0.00 

0.00 

108.34 

528.22) 

2674774.29] 

0.00 

0.00 

120.10 

642.42) 

2665280.29] 

0.00 

0.00 

132.04 

768.47) 

2655786.29] 

0.00 

0.00 

144.18 

906.571 

2646292.29] 

0.00 

0.00 

156.50 

1056.89? 

2636798.29] 

0.00 

0.00 

169.02 

1219.64) 

2627304.29 

0.00 

0.00 

18.1.74 

1395.00 

2617810.29] 

0.00 

0.00 

194.66 

1583.18? 

2608316.29] 

0.00 

0.00 

207.79 

1784.39? 

2598822.29] 

0.00 

0.00 

221.13 

1998.83 

2589328.29] 

0.00 

0.00 

234.69 

2226.72 

2579834.2^] 

0.20 

0.07 

248.46 

2468.28 

2570340.29 

0.81 

0.54 

262.45 

2723.71 

2560846.29 

1.82 

1.81 

276.67 

2993.25 

2551352.29 

3.24 

4.31 

291.10 

3277.12 

2541858.29 

5.08 

8.43 

305.76 

3575.53 

2532364.291 

7.34 

14.61 

320.64 

3888.721 

2522870.29 

10.02 

23.26 

335.74 

4216.88 

2513376.29 

13.14 

34.80 

351.06 

4560.27 

2503882.29 

16.68 

49.67 

366.61 

4919.08 

2494388.29 

20.66 

68.31 

382.37 

5293.55 

2484894.29 

25.09 

91.15 

398.36 

5683.90 

2475400.29 

29.96 

118.64 

414.58 

6090.35 

2465906.29 

35.28 

151.22 

431.01 

6513.13 

2456412.29 

41.06 

189.35 

447.66 

6952.45 

2446918.29 

47.30 

233.49 

464.54 

7408.53 

2437424.29 

54.00 

284.10 

481.64 

7881.60 

2427930.29 

61.17 

341.64 

498.95 

8371.88 

2418436.29 

68.82 

406.6C 

516.49 

8879.58 

2408942.29 

76.94 

479.44 

534.24 

9404.93 

2399448.29 

85.56 

560.6! 

552.21 

9948.13 

2389954.29 

94.66 

650.71 570.39 

10509.42 

2380460.29 

104.2! 

750.12 

> 588.79 

11088.99 

2370966.29 

114.3! 

859.36 

3 607.4C 

11687.0' 

2361472.29 

124.9^ 

1 978.96 

3 626.22 

12303.8) 

5 2351978.29 

136.0! 

j 1109.4* 

t 645.2! 

5 12939.5) 

J 2342484.29 

147.6) 

3 1251.21 

S 664.4) 

5 13594.42 

> 2332990.29 

159.8: 

2 1404.9' 

7 683.91 

L 14268.6) 

5 2323496.29 

172.4< 

) 1571.0' 

7 703.5! 

j 14962.32 

> 2314002.29 

185.6) 

3 1750.1: 

2 723.3! 

) 15675.7' 

7 2304508.29 

199.4: 

2 1942.6: 

2 743.42 

2 16409.11 

S 2295014.29 

213.6! 

9 2149.1! 

3 763.6- 

i\ 17162.6' 

7 2285520.29 



l/AE448data 






58 










68 



7 8 



82 


3 




8 6 


87 


8 8 


8 9 




09 


92 








98 






52.00 


53.00 


54.00 


55.00 


56.00 


57.00 


58.00 


59.00 


60.00 


61.00 


62.00 


63.00 


64.00 


65.00 


66.00 


67.00 


68.00 


69.00 


70.00 


71.00 


72.00 


73.00 


74.00 


75.00 


76.00 


77.00 


78.00 


79.00 


80.00 


81.00 


82.00 


83.00 


84.00 


85.00 


86.00 


87.00 


88.00 


89.00 


90.00 


91.00 


92.00 


93.00 


94.00 


95.00 


96.00 


97.00 


98.00 


99.00 


100.00 


101.00 


102.00 


103.00 


104.00 


B __ 


228.50 


243.87 


259.79 


276.26 


293.31 


310.92 


329.10 


347.87 


367.21 


387.15 


407.67 


428.79 


450.52 


472.84 


495.78 


519.33 


543.50 


568.29 


593.71 


619.75 


646.57 


674.10 


702.30 


731.19 


760.77 


791.04 


822.01 


853.68 


886.05 


919.12 


952.91 


987.41 


1022.63 


1058.57 


1095.14 


1132.43 


1170.51 


1209.39 


1249.06 


1289.53 


1330.80 


1372.87 


1415.74 


1459.41 


1503.90 


1549.18 


1595.28 


1642.19 


1689.91 


1738.44 


1787.78 


1837.94 


1888.92 


C _ 


2370.17 


2606.31 


2858.09 


3126.07 


3410.81 


3712.87 


4032.84 


4371.28 


4728.77 


5105.90 


5503.26 


5921.44 


6361.05 


6822.68 


7306.94 


7814.45 


8345.81 


8901.66 


9482.60 


10089.28 


10722.38 


11382.66 


12070.80 


12787.49 


13533.41 


14309.26 


15115.72 


15953.51 


16823.31 


17725.83 


18661.79 


19631.90 


20636.86 


21677.40 


22754.22 


23867.94 


25019.34 


26209.23 


27438.39 


28707.62 


30017.71 


31369.48 


32763.71 


34201.22 


35682.81 


37209.28 


38781.45 


40400.11 


42066.09 


43780.20 


45543.24 


47356.03 


49219.40 



784.05 


804.65 


825.42 


846.38 


867.51 


888.81 


910.28 


931.90 


953.69 


975.62 


997.71 


1019.93 


1042.29 


1064.78 


1087.40 


1110.14 


1132.99 


1155.95 


1179.02 


1202 .18 


1225.72 


1249.47 


1273.38 


1297.42 


1321.60 


1345.89 


1370.30 


1394.82 


1419.43 


1444.12 


1468.90 


1493.74 


1518.64 


1543.59 


1568.45 


1593.33 


1618.33 


1643.44 


1668.64 


1693.93 


1719.28 


1744.69 


1770.14 


1795.63 


1821.12 


1846.63 


1872.13 


1897.60 


1923.05 


1948.46 


1973.81 


1999.10 


2024.31 


S _ 


17936.50 


18730.83 


19545.85 


20381.74 


21238.67 


22116.82 


23016.35 


23937.43 


24880.21 


25844.85 


26831.51 


27840.31 


28871.41 


29924.94 


31001.02 


32099.78 


33221.33 


34365.80 


35533.27 


36723.86 


37937.77 


39175.35 


40436.77 


41722.16 


43031.66 


44365.40 


45723.49 


47106.04 


48513.15 


49944.92 


51401.43 


52882.74 


54388.92 


55920.03 


57476.08 


59056.96 


60662.78 


62293.65 


63949.69 


65630.97 


67337.57 


69069.56 


70826.97 


72609.86 


74418.23 


76252.11 


78111.49 


79996.35 


81906.68 


83842.44 


85803.58 


87790.04 


898C1.75 


2276026.29 


2266532.29 


2257038.29 


2247544.29 


2238050.29 


2228556.29 


2219062.29 


2209568.29 


2200074.29 


2190580.29 


2181086.29 


2171592.29 


2162098.29 


2152604.29 


2143110.29 


2133616.29 


2124122.29 


2114628.29 


2105134.29 


2095640.29 


2086146.29 


2076652.29 


2067158.29 


2057664.29 


2048170.29 


2038676.29 


2029182.29 


2019688.29 


2010194.29 


2000700.29 


1991206.29 


1981712.29 


1972218.29 


1962724.29 


1953230.29 


1943736.29 


1934242.29 


1924748.29 


1915254.29 


1905760.29 


1896266.29 


1886772.29 


1877278.29 


1867784.29 


1858290.29 


1848796.29 


1839302.29 


1829808.29 


1820314.29 


1810820.29 


1801326.29 


1791832.29 


1782338.29 
























l/AE 448 data 




L 07 


8 


9 





112 







118 











128 








138 










148 


9 










158 



105.00 


106.00 


. 107.00 


108.00 


109.00 


110.00 


111.00 


112.00 


113.00 


114.00 


115.00 


116.00 


117.00 


118.00 


119.00 


120.00 


121.00 


122.00 


123.00 


124.00 


125.00 


126.00 


127.00 


128.00 


129.00 


130.00 


131.00 


132.00 


133.00 


134.00 


135.00 


136.00 


137.00 


138.00 


139.00 


140.00 


141.00 


142.00 


143.00 


144.00 


145.00 


146.00 


147.00 


148.00 


149.00 


150.00 


151.00 


152.00 


153.00 


154.00 


155.00 


156.00 


157.00 


B 


1940.72 


1993.34 


2046.78 


2101.05 


2156.14 


2212.06 


2268.81 


2326.39 


2384.81 


2444.07 


2504.16 


2565.10 


2626.88 


2689.52 


2753.00 


2817.33 


2882.53 


2948.58 


3015.50 


3083.28 


3151.94 


3221.47 


3291.88 


3363.17 


3435.34 


3508.41 


3582.37 


3657.23 


3733.00 


3809.67 


3887.26 


3965.76 


4045.19 


4125.54 


4206.83 


4289.05 


4372.01 


4455.50 


4539.54 


4624.12 


4709.25 


4794.94 


4881.21 


4968.05 


5007.92 


5047.90 


5088.01 


5128.23 


5168.57 


5209.03 


5249.62 


5290.32 


5331.15 


C 


51134.15 


53101.11 


55121.10 


57194.95 


59323.47 


61507.50 


63747.86 


66045.40 


68400.93 


70815.30 


73289.34 


75823.90 


78419.83 


81077.96 


83799.14 


86584.24 


89434.10 


92349.58 


95331.55 


98380.86 


101498.40 


104685.03 


107941.63 


111269.08 


114668.26 


118140.06 


121685.38 


125305.10 


129000.14 


132771.40 


136619.79 


140546.22 


144551.61 


148636.90 


152803.00 


157050.86 


161381.35 


165795.06 


170292.53 


174874.31 


179540.95 


184293.00 


189131.03 


194055.61 


199043.58 


204071.48 


209139.43 


214247.54 


219395.93 


224584.72 


229814.04 


235084.00 


240394.72 


2049.44 


2074.47 


2099.40 


2124.20 


2148.89 


2173.43 


2197.83 


2222.07 


2246.15 


2270.06 


2293.78 


2317.30 


2340.63 


2363.74 


2386.64 


2409.31 


2431.74 


2453.93 


2475.86 


2497.53 


2518.93 


2540.05 


2560.88 


2581.42 


2601.65 


2621.57 


2641.16 


2660.43 


2679.35 


2697.93 


2716.16 


2734.02 


2751.50 


2768.61 


2785.32 


2801.63 


2817.90 


2834.48 


2851.38 


2868.61 


2886.17 


2904.05 


2922.28 


2940.84 


2932.30 


2923.83 


2915.44 


2907.13 


2898.90 


2890.74 


2882.66 


2874.66 


2866.74 


91838.63 


93900.60 


95987.54 


98099.35 


100235.90 


102397.07 


104582.72 


106792.68 


109026.81 


111284.93 


113566.86 


115872.41 


118201.40 


120553.60 


122928.81 


125326.81 


127747.35 


130190.21 


132655.12 


135141.84 


137650.09 


140179.60 


142730.09 


145301.26 


147892.82 


150504.45 


153135.84 


155786.66 


158456.58 


161145.26 


163852.33 


166577.45 


169320.24 


172080.33 


174857.33 


177650.84 


180460.57 


183286.74 


186129.64 


188989.61 


191866.97 


194762.05 


197675.19 


200606.72 


203543.29 


206471.35 


209390.98 


212302.26 


215205.27 


218100.08 


220986.77 


223865.42 


226736.11 


1772844.29 


1763350.29 


1753856.29 


1744362.29 


1734868.29 


1725374.29 


1715880.29 


1706386.29 


1696892.29 


1687398.29 


1677904.29 


1668410.29 


1658916.29 


1649422.29 


1639928.29 


1630434.29 


1620940.29 


1611446.29 


1601952.29 


1592458.29 


1582964.29 


1573470.29 


1563976.29 


1554482.29 


1544988.29 


1535494.29 


1526000.29 


1516506.29 


1507012.29 


1497518.29 


1488024.29 


1478530.29 


1469036.29 


1459542.29 


1450048.29 


1440554.29 


1431060.29 


1421566.29 


1412072.29 


1402578.29 


1393084.29 


1383590.29 


1374096.29 


1364602.29 


1190804.69 


1187007.09 


1183209.49 


1179411.89 


1175614.29 


1171816.69 


1168019.09 


1164221.49 


1160423.89 


























































l / AE 448 data 



L 6 0 



158.00 


159.00 


160.00 


161.00 


162.00 


163.00 


164.00 


165.00 


166.00 


167.00 


168.00 


169.00 


170.00 


171.00 


172.00 


173.00 


174.00 


175.00 


178 

176.00 

ESH 

177.00 

180 

178.00 

181 

179.00 

182 

180.00 

183 

181.00 

184 

182.00 

185 

183.00 

186 

184.00 

187 

185.00 

188 

186.00 

189 

187.00 

Em 

roi 

188.00 

189.00 



190.00 


191.00 


192.00 


193.00 


194.00 


195.00 


196.00 


197.00 


198.00 


199.00 


200.00 


201.00 


202.00 


203.00 


204.00 


205.00 


206.00 


207.00 


208.00 


209.00 


210.00 


B _ 


5372.10 


5413.17 


5454.37 


5495.79 


5537.55 


5579.64 


5622.06 


5664.81 


5707.89 


5751.29 


5795.02 


5839.08 


5883.45 


5928.15 


5973.17 


6018.51 


6064.17 


6110.14 


6156.43 


6203.04 


6249.95 


6297.18 


6344.72 


6392.57 


6440.72 


6489.18 


6537.94 


6587.01 


6636.37 


6686.04 


6736.00 


6786.26 


6836.81 


6887.66 


6938.80 


6990.23 


7041.94 


7093.94 


7146.23 


7198.79 


7251.64 


7304.72 


7357.99 


7411.46 


7465.11 


7518.97 


7573.01 


7627.26 


7681.71 


7736.36 


7791.21 


7846.27 


7901.53 


245746.33 


251138.95 


256572.71 


262047.76 


267564.41 


273122.98 


278723.80 


284367.21 


290053.53 


295783.09 


301556.22 


307373.25 


313234.48 


319140.26 


325090.89 


331086.71 


337128.02 


343215.15 


349348.41 


355528.12 


361754.59 


368028.13 


374349.06 


380717.67 


387134.29 


393599.21 


400112.75 


406675.19 


413286.86 


419948.04 


426659.03 


433420.14 


440231.65 


447093.87 


454007.07 


460971.56 


467987.62 


475055.53 


482175.59 


489348.08 


496573.27 


503851.44 


511182.78 


518567.49 


526005.76 


533497.78 


541043.75 


548643.88 


556298.34 


564007.36 


571771.13 


579589.85 


587463.73 


2858.89 


2851.13 


2843.45 


2835.67 


2827.60 


2819.25 


2810.61 


2801.68 


2792.45 


2782.92 


2773.08 


2762.94 


2752.48 


2741.71 


2730.63 


2719.22 


2707.48 


2695.42 


2683.02 


2670.28 


2657.21 


2643.79 


2630.02 


2615.90 


2601.43 


2586.60 


2571.41 


2555.85 


2539.92 


2523.61 


2506.93 


2489.87 


2472.42 


2454.59 


2436.36 


2417.74 


2398.72 


2379.29 


2359.46 


2339.21 


2318.55 


2297.70 


2276.90 


2256.14 


2235.42 


2214.75 


2194.13 


2173.55 


2153.01 


2132.52 


2112.07 


2091.68 


2071.32 


229598.92 


232453.93 


235301.21 


238140.79 


240972.45 


243795.90 


246610.86 


249417.03 


252214.12 


255001.83 


257779.85 


260547.89 


263305.62 


266052.75 


268788.95 


271513.89 


274227.27 


276928.75 


279617.99 


282294.67 


284958.44 


287608.97 


290245.91 


292868.90 


295477.60 


298071.64 


300650.68 


303214.34 


305762.25 


308294.04 


310809.35 


313307.78 


315788.96 


318252.50 


320698.01 


323125.10 


325533.36 


327922.40 


330291.80 


332641.17 


334970.09 


337278.21 


339565.51 


341832.03 


344077.80 


346302.89 


348507.33 


350691.16 


352854.43 


354997.19 


357119.49 


359221.36 


361302.85 


1156626.29 


1152828.69 


1149031.09 


1145233.49 


1141435.89 


1137638.29 


1133840.69 


1130043.09 


1126245.49 


1122447.89 


1118650.29 


1114852.69 


1111055.09 


1107257.49 


1103459.89 


1099662.29 


1095864.69 


1092067.09 


1088269.49 


1084471.89 


1080674.29 


1076876.69 


1073079.09 


1069281.49 


1065483.89 


1061686.29 


1057888.69 


1054091.09 


1050293.49 


1046495.89 


1042698.29 


1038900.69 


1035103.09 


1031305.49 


1027507.89 


1023710.29 


1019912.69 


1016115.09 


1012317.49 


1008519.89 


1004722.29 


1000924.69 


997127.09 


993329.49 


989531.89 


985734.29 


981936.69 


978139.09 


974341.49 


970543.89 


966746.29 


962948.69 


959151.09 


25 































l / AE 440 data 


211.00 
212.00 
213.00 
214.00 
215.00 
216.00 
217.00 
218.00 
219.00 
220.00 
221.00 
222.00 
223.00 
224.00 
225.00 
226.00 
227.00 
228.00 
229.00 
230.00 
231.00 
232,00 
233.00 
234.00 
235.00 
236.00 
237.00 
238.00 
239.00 
240.00 
241.00 
242.00 
243.00 
244.00 
245.00 
246.00 
247.00 
248.00 
249.00 
250.00 
251.00 
252.00 
253.00 
254.00 
255.00 
256.00 
257.00 
258.00 
259.00 
260.00 
261.00 


7957.00 

8012.68 

8068.58 

8124.69 
8181.01 

8237.55 
8294.31 

8351.29 
8408.50 

8465.93 

8523.59 
8581.47 

8639.59 

8697.94 

8756.53 
8815.35 
8874.41 
8933.72 

8993.27 
9053.06 

9113.10 
9173.40 

9233.94 
9294.74 

9355.80 

9417.12 

9478.70 

9540.54 
9602.66 
9665.04 

9727.69 

9790.62 
9853.83 
9917.31 
9981.08 

10045.14 

10109.48 
10174.11 
10239.04 

10304.27 
10369.79 

10435.62 
10501.75 
10568.19 

10634.94 
10702.01 
10769.39 

10837.10 
10905.13 

10973.49 
11042.18 


595392.98 

603377.80 
611418.41 
619515.03 
627667.86 

635877.12 
644143.03 

652465.81 

660845.69 

669282.88 

677777.62 

686330.13 
694940.64 
703609.39 
712336.60 


D 

2051.02 

2030.76 

2010.54 

1990.38 

1970.26 

1950.19 
1930.17 

1910.20 
1890.28 

1870.40 
1850.58 
1830.81 
1811.08 

1791.41 
1771.79 
1752.22 


363364.02 

365404.90 

367425.55 
369426.01 
371406.32 

373366.55 
375306.73 

377226.91 
379127.15 

381007.49 
382867.97 
384708.66 
386529.60 
388330.85 

390112.44 

391874.44 


955353.49 

951555.89 

947758.29 

943960.69 
940163.09 

936365.49 

932567.89 

928770.29 

924972.69 
921175.09 

917377.49 

913579.89 

909782.29 

905984.69 
902187.09 

898389.49 


262.00 

263.00 


11111.20 


729967.38 

1732.70 

393616.89 

894591.89 

738871.42 

1713.23 

395339.85 

890794.29 

747834.89 

1693.81 

397043.36 

886996.69 

756858.03 

1674.45 

398727.49 

883199.09 

765941.09 

1655.14 

400392.28 

879401.49 

775084.32 

1635.09 

402037.79 

875603.89 

784287.97 

1616.68 

403664.07 

871806.29 

793552.29 

1597.54 

405271.18 

868008.69 

802877.54 

1578.44 

406859.16 

864211.09 

812263.98 

1559.41 

408428.08 

860413.49 

821711.86 

1540.42 

409977.99 

056615.89 

831221.46 

1521.50 

411508.95 

852818.29 

840793.04 

1502.63 

413021.01 

849020.69 

850426.86 

1483.81 

414514.22 

845223.09 

860123.21 

1465.06 

415988.65 

841425.49 

869882.34 

1446.36 

417444.36 

837627.89 

879704.54 

1427.72 

418881.39 

833830.29 

889590.09 

1409.14 

420299.81 

830032.69 

899539.26 

1390.61 

421699.68 

826235.09 

909552.35 

1372.15 

423081.06 

822437.49 

919629.641 

1353 . 7 S 1 

424444.00 

818639.89 

929771.41 

1335.40 

425788.57 

814842.29 

939977.96 

1317.12 

427114.83 

811044.69 

950249.59 

1298.90 

428422.841 

807247.09 

960586.59 

1280.74 

429712.66 

803449.49 

970989.27 

1262.65 

430984.35 

799651.89 

981457.93 

1244.62 

432237.98 

795854.29 

991992.87 

1226.65 

433473.60 

792056.69 

1002594.41 

1208.74 

434691.29 

788259.09 

1013262.86 

1190.90 

435891.11 

784461.49 

1023998.53 

1173.13 

437073.12 

780663.89 

1034801.75 

1155.42 

438237.38 

776866.29 

1045672.84 

1137.77 

439383.97 

773068.69 

1056612.12 

1120.20 

440512.95 

769271.09 

1067619.92 

1102.69 

441624.39 

765473.49 

1078696.58 

1085.25 

442718.35 

761675.89 

1089842.42 

1067.88 

443794.91 

757878.29 




























l / AE 448 data 




268 



278 



28 


281 


282 



264.00 


265.00 


266.00 


267.00 


268.00 


269.00 


270.00 


271.00 


272.00 


273.00 


274.00 


275.00 


276.00 


277.00 


278.00 


279.00 


280.00 


281.00 


284 

282.00 

285 

283.00 

286 

284.00 

287 

285.00 

288 

286.00 

289 

287.00 


288.00 

fog 

289.00 



318 


290.00 


291.00 


292.00 


293.00 


294.00 


295.00 


296.00 


297.00 


298.00 


299.00 


300.00 


301.00 


302.00 


303.00 


304.00 


305.00 


306.00 


307.00 


308.00 


309.00 


310.00 


311.00 


312.00 


313.00 


314.00 


315.00 


316.00 


B 


11250.27 


11320.31 


11390.71 


11461.46 


11532.57 


11604.03 


11675.86 


11748.06 


11820.62 


11893.57 


11966 . 8 ! 


12040.60 


12114 . 6 ! 


12189 . 1 ! 


12264.06 


12339.35 


12415.04 


12491.14 


12567.66 


12644.59 


12721.95 


12799.74 


12877.96 


12956.62 


13035.73 


13115.29 


13195.30 


13275.77 


13356.70 


13438.11 


13520.00 


13602.37 


13685.22 


13768.57 


13852.43 


13936.79 


14021.66 


14107.05 


14192.97 


14279.42 


14366.42 


14453.96 


14542.05 


14630.70 


14719.93 


14809.72 


14900.11 


14991.08 


15082.65 


15174.84 


15267.63 


15361.05 


15455.11 


1101057.82 


1112343.08 


1123698.56 


1135124.62 


1146621.60 


1158189.87 


1169829.78 


1181541.71 


1193326.02 


1205183.08 


1217113.28 


1229116.99 


1241194.60 


1253346.50 


1265573.09 


1277874.76 


1290251.92 


1302704.98 


1315234.34 


I 1327840.43 


1340523.66 


1353284.47 


1366123.29 


1379040.54 


1392036.68 


1405112.15 


1418267.40 


1431502.90 


1444819.09 


1458216.46 


1471695.48 


1485256.62 


1498900.37 


1512627.23 


1526437.69 


1540332.25 


1554311.43 


1568375.74 


1582525.71 


1596761.86 


2 | 1611084.74 


1625494.88 


1639992.83 


1654579.16 


1669254.43 


1684019.21 


1698874.07 


1713819.62 


1728856.44 


1743985.13 


1759206.31 


1774520.60 


1789928.63 


1050.58 


1033.34 


1016.18 


999.09 


982.08 


965.13 


948.26 


931.46 


914.74 


898.09 


881.52 


865.03 


848.62 


832.28 


816.02 


799.85 


783.75 


767.74 


751.81 


735.96 


720.20 


704.52 


688.93 


673.43 


658.01 


642.69 


627.45 


612.31 


597.26 


582.30 


567.44 


552.67 


537.99 


523.42 


508.95 


494.57 


480.30 


466.13 


452.06 


438.09 


424.24 


410.49 


396.85 


383.32 


369.90 


356.59 


343.40 


330.33 


317.37 


304.53 


291.81 


279.22 


266.75 


444854.13 


445896.09 


446920.84 


447928.48 


448919.06 


449892.65 


450849.34 


451789.20 


452712.30 


453618.71 


454508.51 


455381.78 


456238.60 


457079.04 


457903.18 


458711.11 


459502.91 


460278.64 


461038.41 


461782.28 


462510.36 


463222.71 


463919.43 


464600.60 


465266.31 


465916.66 


466551.72 


467171.59 


467776.37 


468366.14 


468941.00 


469501.04 


470046.36 


470577.06 


471093.24 


471594.99 


472082.41 


472555.61 


473014.70 


473459.76 


473890.92 


474308.27 


474711.93 


475102.01 


475478.61 


475841.85 


476191.83 


476528.69 


476852.53 


477163.47 


477461.64 


477747.14 


478020.11 


754080.69 


750283.09 


746485.49 


742687.89 


738890.29 


735092.69 


731295.09 


727497.49 


723699.89 


719902.29 


716104.69 


712307.09 


708509.49 


704711.89 


700914.29 


697116.69 


693319.09 


689521.49 


685723.89 


681926.29 


678128.69 


674331.09 


670533.49 


666735.89 


662938.29 


659140.69 


655343 . 0 ! 


651545.49 


647747.89 


643950.29 


640152.69 


636355.09 


632557.49 


628759 . 8 ! 


624962.29 


621164.69 


617367.09 


613569.49 


609771.89 


605974.29 


602176.69 


598379.09 


594581.49 


590783.89 


586986.29 


583188.69 


579391.09 


575593.49 


571795.89 


567998.29 


564200.69 


560403.09 


556605.49 
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l / AE 448 data 



32 


322 



328 





332 




335 




348 


9 


35 


35 


352 


353 


355 


356 


358 


36 


362 


363 


365 


368 


317.00 


318.00 


319.00 


320.00 


321.00 


322.00 


323.00 


324.00 


325.00 


326.00 


327.00 


328.00 


329.00 


330.00 


331.00 


332.00 


333.00 


334.00 


335.00 


338 

336.00 

339 

337.00 

EED 

BED 

338.00 

339.00 


340.00 


341.00 


342.00 


343.00 


344.00 


345.00 


346.00 


347.00 


348.00 


349.00 


350.00 


351.00 


352.00 


353.00 


354.00 


355.00 


356.00 


357.00 


358.00 


359.00 


360.00 


361.00 


362.00 


363.00 


364.00 


365.00 


366.00 


367.00 


368.00 


369.00 


B 


15549.80 


15645.15 


15741.15 


15837.83 


15935.18 


16033.22 


16131.96 


16231.41 


16331.57 


16432.47 


16534.10 


16636.49 


16739.64 


16843.56 


16948.27 


17053.78 


17160.10 


17267.24 


17375.21 


17484.04 


17593.73 


17704.29 


17815.75 


17928.11 


18041.39 


18155.61 


18270.78 


18386.91 


18504.03 


18622.15 


18741.29 


18861.47 


18982.70 


19105.00 


19228.39 


19352.90 


19478.54 


19605.33 


19733.29 


19862.46 


19992.84 


20124.46 


20257.35 


20391.53 


20527.03 


20663.87 


20802.07 


20941.68 


21082.71 


21225.20 


21369.17 


21514.65 


21661.69 


1805431.03 


1821028.45 


1836721.54 


1852510.98 


1868397.42 


1884381.57 


1900464.10 


1916645.72 


1932927.15 


1949309.11 


1965792.33 


1982377.57 


1999065.57 


2015857.11 


2032752.96 


2049753.92 


2066860.79 


2084074.38 


2101395.54 


2118825.09 


2136363.90 


2154012.84 


2171772.78 


2189644.63 


2207629.31 


2225727.73 


2243940.84 


2262269.60 


2280714.99 


I 2299278.00 


2317959.63 


2336760.93 


2355682.92 


2374726.60 


2393893.28 


2413183.83 


2432599.45 


2452141.29 


2471810.50 


2491608.27 


2511535.82 


2531594.36 


2551785.16 


2572109.49 


2592568.66 


2613163.99 


2633896.84 


2654768.60 


2675780.68 


2696934.51 


2718231.56 


2739673.35 


2761261.39 


254.40 


242.18 


230.09 


218.13 


206.30 


194.61 


183.05 


171.64 


160.36 


149.22 


138.23 


127.38 


116.69 


106.14 


95.75 


85.51 


75.43 


65.51 


55.75 


46.16 


36.74 


27.49 


18.41 


9.50 


0.78 


- 7.77 


- 16.12 


- 24.29 


- 32.27 


- 40.06 


- 47.65 


- 55.03 


- 62.21 


- 69.18 


- 75.94 


- 82.49 


- 88.81 


- 94.91 


- 100.79 


- 106.43 


- 111.83 


- 116.99 


- 121.91 


- 126.58 


- 130.99 


- 135.14 


- 139.02 


- 142.64 


- 145.97 


- 149.03 


- 151.79 


- 154.27 


- 156.44 


478280.68 


478528.96 


478765.08 


478989.18 


479201.39 


479401.83 


479590.65 


479767.98 


479933.97 


480088.74 


480232.46 


480365.25 


480487.28 


480598.68 


480699.61 


480790.23 


480870.68 


480941.14 


481001.76 


481052.71 


481094.14 


481126.24 


481149.17 


481163.11 


481168.24 


481164.73 


481152.77 


481132.54 


481104.24 


481068.06 


481024.19 


480972.83 


480914.20 


480848.48 


480775.90 


480696.66 


480610.99 


480519.11 


480421.24 


480317.61 


480208.46 


480094.03 


479974.56 


479850.29 


479721.49 


479588.40 


479451.30 


479310.45 


479166.12 


479018.59 


478868.16 


478715.10 


478559.73 


552807.89 


549010.29 


545212.69 


541415.09 


537617.49 


533819.89 


530022.29 


526224.69 


522427.09 


518629.49 


514831.89 


511034.29 


507236.69 


503439.09 


499641.49 


495843.89 


492046.29 


488248.69 


484451.09 


480653.49 


476855.89 


473058.29 


469260.69 


465463.09 


461665.49 


457867.89 


454070.29 


450272.69 


446475.09 


442677.49 


438879.89 


435082.29 


431284.69 


427487.09 


423689.49 


419891.89 


416094.29 


412296.69 


408499.09 


404701.49 


400903.89 


397106.29 


393308.69 


389511.09 


385713.49 


381915.89 


378118.29 


374320.69 


370523.09 


366725.49 


362927.89 


359130.29 


355332.69 
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l/AE448data 


1 


1 

A 1 

B [ 

C 

D 

B 

F | 

3 72 

370.00 

21810.31 

2782997.25 

- 158.30 

478402.33 

351535 . 09 ] 

iTtT 

371.00 

21960.54 

2804882.54 

- 159.85 

478243.23 

347737 . 49 ] 

ran 

372.00 

22112.43 

2826918.89 

- 161.08 

478082.74 

343939.89 

CMMlm 

ran 

373.00 

22266.00 

2849107.96 

- 161.98 

477921.18 

340142.29 

nn 

ran 

374.00 

22421.31 

2871451.47 

- 162.55 

477758.89 

336344 . 69 ] 

■ ■! a 

ran 

375.00 

22578.37 

2893951.17 

- 162.77 

477596.20 

332547.09 


376.00 

22737.25 

2916608.82 

- 162.65 

477433.46 

328749.49 

ran 

377.00 

22897.97 

2939426.28 

- 162.16 

477271.02 

324951.89 

■ mJ 

378.00 

23060.59 

2962405.40 

- 161.31 

477109.26 

321154 . 29 ] 

3 81 

379.00 

23225.14 

2985548.10 

- 160.08 

476948.54 

317356 . 69 ] 

3 82 

380.00 

23391.67 

3008856.33 

- 158.46 

476789.23 

313559.09 

383 

381.00 

23560.24 

3032332.12 

- 156.46 

476631.74 

309761 . 49 ] 

3 8 4 

382.00 

23730.89 

3055977.50 

- 154.04 

476476.45 

305963.89 

3 85 

383.00 

23903.67 

3079794.60 

- 151.22 

476323.79 

302166 . 29 ] 

386 

384.00 

24078.64 

3103785.57 

- 147.96 

476174.16 

298368 . 69 ] 

3 87 

385.00 

24255.85 

3127952.63 

- 144.28 

476028.01 

294571 . 09 ] 

3 8 8 

386.00 

24435.37 

3152298.05 

- 140.14 

475885.76 

290773.49 

3 8 9 

387.00 

24617.25 

3176824.16 

- 135.55 

475747.87 

286975.89 

ran 

388.00 

24801.55 

3201533.35 

- 130.48 

475614.82 

283178.29 

iclrll 

389.00 

24988.35 

3226428.09 

- 124.94 

475487.07 

279380 . 69 ] 


390.00 

25177.70 

3251510.89 

- 118.89 

475365.11 

275583.09 


391.00 

25369.68 

3276784.36 

- 112.33 

475249.45 

271785 . 49 ] 

ran 

392.00 

25564.37 

3302251.16 

- 105.25 

475140.62 

267987.89 


393.00 

25761.83 

3327914.02 

- 97.63 

475039.13 

264190.29 

U 9 6 

394.00 

25962.16 

3353775.78 

- 89.46 

474945.54 

260392.69 

km 

395.00 

25962.16 

3379737.94 

- 120.22 

474840.70 

260392 . 69 ] 

3 98 

396.00 

25962.16 

3405700.10 

- 150.98 

474705.10 

260392 . 69 ] 

399 

397.00 

25962.16 

3431662.26 

- 181.74 

474538.74 

260392 . 69 ] 

Inn 

398.00 

25962.16 

3457624.42 

- 212.51 

474341.61 

260392.69 

lEEH 

399.00 

25962.16 

3483586.58 

- 243.27 

474113.72 

260392.691 

Iran 

400.00 

25962.16 

3509548.74 

- 274.04 

473855.07 

260392.69 

Iran 

401.00 

25962.16 

3535510.91 

- 304.80 

473565.64 

260392 . 69 ] 

Inn 

402.00 

25962.16 

3561473.07 

- 335.57 

473245.46 

260392.69 

|nn 

403.00 

25962.16 

3587435.23 

- 366.34 

472894.50 

260392.69 

Inn 

404.00 

25962.16 

3613397.39 

- 397.11 

472512.78 

260392 . 6 ^] 

Iran 

405.00 

25962.16 

3639359.55 

- 427.88 

472100.28 

260392.69 


406.00 

25962.16 

3665321.71 

- 458.65 

471657.02 

260392.69 

IEEE 

407.00 

25962.16 

3691283.87 - 489.42 

471182.98 

260392.69 

Imn 

408.00 

25962.16 

3717246.03 

- 520.20 

470678.17 

260392.69 

Inn 

409.00 

25962.16 

3743208.19 

- 550.97 

470142.58 

260392.69 

IEEE 

410.00 

25962.16 

3769170.36 

- 581.75 

469576.22 

260392.69 

Imn 

411.00 

25962.16 

3795132.52 

- 612.53 

468979.08 

260392.69 

IEEE] 

412.00 

25962.16 

3821094.68 

- 643.31 

468351.16 

260392.69 

IEEE 

413.00 

25962.16 

3847056.84 

- 674.09 

467692.46 

260392.69 

Imn 

414.00 

25962.16 

3873019 . 0 C 

- 704.88 

467002.98 

260392.69 

Inn 

415.00 

25962.16 

3898981.16 

- 735.66 

466282.71 

260392.69 


416.00 

25962.16 

3924943.32 

- 766.45 

465531.65 

260392.69 

Inn 

417.00 

25962.16 

3950905.48 

- 797.24 

464749.80 

260392.69 

Imn 

418.00 

25962.16 

3976867.64 

- 828.04 

463937.16 

260392.69 

Iren 

419 . 0 C 

25962.16 

4002829.81 

- 858.83 

463093.73 

260392 . 69 ] 


420 . 0 C 

25962.16 

4028791 . 9 ' 

- 889.63 

462219.50 

260392.691 
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18 





23 




2 6 


Inc 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


8 

90.00 

9 

90.00 


90.00 




52 


53 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


89.50 


89.00 


88.50 


88.00 


87.50 


87.00 


86.50 


86.00 


85.50 


85.00 


84.50 


84.00 


83.50 


83.00 


82.50 


82.00 


81.50 


81.00 


80.50 


80.00 


79.50 


79.00 


78.50 


78.00 


77.50 


77.00 


76.50 


76.00 


75.50 


75.00 


74.50 


H 


Vtotal 


0.00 


10.13 


20.40 


30.83 


41.41 


52.15 


63.05 


74.11 


85.35 


96.75 


108.34 


120.10 


132.04 


144.18 


156.50 


169.02 


181.74 


194.66 


207.79 


221.13 


234.69 


248.46 


262.45 


276.67 


291.12 


305.80 


320.72 


335.89 


351.31 


366.99 


382.93 


399.15 


415.66 


432.45 


449.54 


466.94 


484.65 


502.69 


521.05 


539.75 


558.80 


578.19 


597.95 


618.07 


638.56 


659.43 


680.69 


702.34 


724.39 


746.84 


769.70 


792.97 


I 


Gamma 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


90.00 


89.95 


89.82 


89.62 


89.36 


89.05 


88.69 


88.29 


87.86 


87.39 


86.91 


86.40 


85.87 


85.32 


84.76 


84.19 


83.60 


83.01 


82.41 


81.80 


81.19 


80.58 


79.96 


79.34 


78.72 


78.09 


77.47 


76.85 


76.22 


75.60 


74.98 


74.37 


Thrust 


0.00 


3635328.12 


3635599.61 


3636068.00 


3636735.59 


3637604.52 


3638676.70 


3639953.82 


3641437.35 


3643128.53 


3645028.31 


3647137.44 


3649456.35 


3651985.22 


3654723.95 


3657672.12 


3660829.04 


3664193.67 


3667764.69 


3671540.44 


3675518.93 


3679697.82 


3684074.43 


3688645.66 


3693408.02 


3698357.60 


3703490.10 


3708800.82 


3714284.68 


3719936.20 


3725749.52 


3731718.41 


3737836.31 


3744096.28 


3750491.08 


3757013.13 


3763654.57 


3770407.26 


3777262.78 


3784212.51 


3791247.58 


3798358.96 


3805537.43 


3812773.65 


3820058.15 


3827381.39 


3834733.78 


3842105.67 


3849487.46 


3856869.55 


3864242.41 


3871596.60 



14.700000 


14.697595 


14.689235 


14.674810 


14.654251 


14.627492 


14.594474 


14.555144 


14.509458 


14.457377 


14.398872 


14.333921 


14.262509 


14.184631 


14.100290 


14.009499 


13.912280 


13.808664 


13.698693 


13.582416 


13.459897 


13.331205 


13.196425 


13.055651 


12.908992 


12.756567 


12.598509 


12.434962 


12.266083 


12.092042 


11.913017 


11.729202 


11.540798 


11.348018 


11.151087 


10.950237 


10.745710 


10.537757 


10.326637 


10.112616 


9.895967 


9.676969 


9.455903 


9.233060 


9.008730 


8.783206 


8.556786 


8.329764 


8.102438 


7.875102 


7.648051 


7.421574 


Den sit: 


0.002377 


0.002376 


0.002375 


0.002373 


0.002369 


0.002365 


0.002360 


0.002353 


0.002346 


0.002338 


0.002328 


0.002318 


0.002306 


0.002293 


0.002280 


0.002265 


0.002249 


0.002233 


0.002215 


0.002196 


0.002176 


0.002155 


0.002134 


0.002111 


0.002087 


0.002063 


0.002037 


0.002011 


0.001983 


0.001955 


0.001926 


0.001896 


0.001866 


0.001835 


0.001803 


0.001771 


0.001737 


0.001704 


0.001670 


0.001635 


0.001600 


0.001565 


0.001529 


0.001493 


0.001457 


0.001420 


0.001384 


0.001347 


0.001310 


0.001273 


0.001237 


0.001200 
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74.00 

73.50 
73.00 

72.50 
72.00 

71.50 
71.00 

70.50 
70.00 

69.50 
69.00 

68.50 

68.00 

67.50 
67.00 

66.50 

66.00 

65.50 
65.00 

64.50 
64.00 

63.50 
63.00 

62.50 
62.00 

61.50 
61.00 

60.50 
60.00 

59.50 
59.00 

58.50 
58.00 

57.50 
57.00 

56.50 
56.00 

55.50 
55.00 

54.50 
54.00 

53.50 
53.00 

52.50 
52.00 

51.50 
51.00 

50.50 
50.00 

49.50 
49.00 

48.50 
48.00 


H 

? 16.67 
840.79 
865.34 
890.33 
915.75 
941.62 
967.94 
994.71 
1021.94 
1049.63 
1077.78 
1106.40 
1135.49 
1165.05 
1195.09 

1225.61 

1256.61 
1288.09 
1320.06 
1352.53 

1385.80 
1419.72 
1454.21 

1489.28 
1524.92 

1561.15 

1597.95 
1635.32 

1673.28 

1711.81 

1750.91 
1790.60 
1830.86 
1871.70 

1912.95 
1954.76 
1997.27 
2040.47 
2084.35 

2128.92 

2174.16 
2220.07 
2266.66 
2313.91 
2361.82 
2410.40 
2459.63 
2509.52 
2560.06 
2611.25 
2663.10 
2715.59 
2768.73 


3878922.82 

3886211.90 
3893454.85 

3900642.90 

3907767.51 
3914820.40 
3921793.57 
3928679.34 
3935470.33 
3942159.56 
3948740.37 

3955206.52 

3961552.16 
3967771.84 
3973860.55 
3979813.71 

3985627.17 
3991297.25 
3996820.67 
4002194.65 
4035115.95 
4045271.09 
4055118.40 

4064655.73 

4073881.73 

4082795.75 
4091397.87 
4099688.91 
4107670.35 
4115344.34 
4122713.67 

4129781.75 
4136552.54 
4143030.57 

4134351.81 
4147546.01 

4160122.37 

4172095.81 

4183481.81 

4194296.37 
4204555.91 
4214277.28 


K 

7.195959 
6.971488 
6.748437 
6.527077 
6.307670 
6.090473 
5.875730 
5.663679 
5.454546 
5.248548 
5.045888 
4.846759 
4.651342 
4.459803 
4.272298 
4.088967 
3.909938 


L 

0.001163 
0.001127 
0.001091 
0.001055 
0.001020 
0,000985 
0,000950 
0,000916 
0.000882 


0.000849 

0.000816 

0,000784 

0,000752 

0.000721 

0,000691 

0,000661 

0.000632 


3.7353251 

0.000604 

3.565228 

0 . 000576 ] 

3.3997331 

0 . 000550 ] 

3.238903 

0 . 000524 ] 

3.082764 

0 . 000498 ] 

2.931357 

0 . 000474 ] 

2.7847161 

0 . 000450 ] 

2.6428631 

0 . 000427 ] 

2.505806 

0 . 000405 ] 

2.3735441 

0 . 000384 ] 

2.246066 

0 . 000363 ] 

2.123348 

0.000343 

2.005357 

0.0003241 

1.892050 

0.0003061 

1.783376 

0.000288 

1.679272 

0.000272 

1.579669 

0 . 000255 ] 

1.484493 

0 . 000240 ] 

1.393666 

0.000225 

1.307092 

0 .00021 l ] 

1.224668 

0 . 000198 ] 

1.146289 

0.000185 

1.071843 

0 . 000173 ] 

1.001217 

0 . 000162 ] 

0.934297 

0 . 00015 l ] 

0.870963 

0.000141 

0.811096 

0 . 00013 l ] 

0.754574 

0 .000122] 

0.701275 

0 . 000113 ] 

0.651076 

0 . 000105 ] 

0.603855 

0 . 00009 ^ 

0.559488 

0.000090 

) 0.5178541 

0 . 000084 ] 

> 0-478832 

0.000077 

3 0.442302 

0.0000 721 

0.4081461 

0 . 000066 ] 



l / AE 448 data 




108 


9 









117 


118 



120 




138 







148 


9 




152 






158 



47.50 


47.00 


46.50 


46.00 


45.50 


45.00 


44.50 


44.00 


43.50 


43.00 


42.50 


42.00 


41.50 


41.00 


40.50 


40.00 


39.50 


39.00 


38.50 


38.00 


37.50 


37.00 


36.50 


36.00 


35.50 


35.00 


34.50 


34.00 


33.50 


33.00 


32.50 


32.00 


31.50 


31.00 


30.50 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


30.00 


H _ 


2822.52 


2876.95 


2932.03 


2987.75 


3044.12 


3101.13 


3158.79 


3217.10 


3276.05 


3335.66 


3395.91 


3456.82 


3518.39 


3580.61 


3643.49 


3707.04 


3771.25 


3836.13 


3901.68 


3967.91 


4034.81 


4102.40 


4170.68 


4239.65 


4309.31 


4379.67 


4450.74 


4522.52 


4595.02 


4668.24 


4742.18 


4816.86 


4892.27 


4968.43 


5045.33 


5123.00 


5201.44 


5280.70 


5360.76 


5441.63 


5523.31 


5605.80 


5689.10 


5773.22 


5803.24 


5833.54 


5864.10 


5894.93 


5926.02 


5957.38 


5989.01 


6020.89 


6053.04 


46.56 


46.14 


45.73 


45.31 


44.90 


44.50 


44.09 


43.69 


43.28 


42.89 


42.49 


42.09 


41.70 


41.31 


40.92 


40.54 


40.15 


39.77 


39.39 


39.01 


38.63 


38.25 


37.88 


37.51 


37.14 


36.77 


36.40 


36.03 


35.67 


35.31 


34.94 


34.58 


34.22 


33.87 


33.51 


33.15 


32.80 


32.46 


32.13 


31.81 


31.50 


31.20 


30.91 


30.62 


30.35 


30.08 


29.81 


29.55 


29.29 


29.03 


28.77 


28.52 


28.27 


4295343.40 


4299665.45 


4303692.05 


4307438.91 


4310921.32 


4314154.13 


4317151.69 


4319927.87 


4322496.04 


4324869.03 


4327059.16 


4329078.22 


4330937.47 


4332647.63 


4334218.91 


4335661.00 


4336983.07 


4338193.79 


4339301.36 


4340313.48 


4341237.42 


4342079.97 


4342847.52 


4343546.04 


4344181.09 


4344757.88 


4345281.22 


4345755.63 


4346185.26 


4346573.98 


4346925.35 


4347242.69 


4347529.02 


4347787.15 


4348019.66 


4348228.91 


4348417.09 


4348586.19 


4348738.07 


4348874.41 


4348996.71 


4349106.36 


4349204.59 


4349292.55 


1739748.39 


1739776.21 


1739800.88 


1739822.77 


1739842.21 


1739859.47 


1739874.80 


1739888.42 


1739900.52 


0.376249 


0.346496 


0.318778 


0.292985 


0.269012 


0.246758 


0.226123 


0-207012 


0.189334 


0.172998 


0.157922 


0.144023 


0.131224 


0.119451 


0.108635 


0.098708 


0.089607 


0.081272 


0.073648 


0.066681 


0.060320 


0.054520 


0.049237 


0.044428 


0.040057 


0.036086 


0.032483 


0.029218 


0.026260 


0.023584 


0.021165 


0.018981 


0.017010 


0.015233 


0.013632 


0.012192 


0.010897 


0.009732 


0.008687 


0.007748 


0.006907 


0.006152 


0.005475 


0.004870 


0.004330 


0.003851 


0.003427 


0.003050 


0.002716 


0.002419 


0.002155 


0.001920 


0.001712 


0.000061 


0.000056 


0.000052 


0.000047 


0.000044 


0.000040 


0.000037 


0.000033 


0.000031 


0.000028 


0.000026 


0.000023 


0.000021 


0.000019 


0.000018 


0.000016 


0.000014 


0.000013 


0.000012 


0.000011 


0.000010 


0.000009 


0.000008 


0.000007 


0.000006 


0.000006 


0.000005 


0.000005 


0.000004 


0.000004 


0.000003 


0.000003 


0.000003 


0.000002 


0.000002 


0.000002 


0.000002 


0.000002 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 
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l / AE 448 data 


160 

1611 

162 

163 p 

164 

i 6 5 r 
16 6 ] 

1 6 7 ~ 

1 6 81 

I 6 9| 

17 0 I 

171~ 

172 

17 31 

17 4 
175) 
17 6 

17 7 

178 

179 

1 8 01 

181 

1 8 21 

183f 

18 4] 

1 8 5 j 

I I 8 61 

1 1 8 7 r 
188 
189 

119 o r 

11 91 j 

192] 

193 

194 
[1 9 5] 

196] 

197 

198 

199 

200 
201 

2 02 
2 03 
204 
2 05 
206 

207 

208 

209 

210 

till 

212 


30.001 

30.00 

30.00 

29.50 
29.00 ~ 

28.50 
28.00 

27.50 
27.00 

26.50 
26.00 

25.50 
25.00 ~ 

24.50 
24.00 

23.50 
23.00 

22.50 

22.00 

21.50 

21.00 

20.50 

20.00 

19.50 
19.00 

18.50 
18.00 

17.50 
17 . 00 [ 

16.50 
16.00 

15.50 
15.00 
14 . 50 ] 
14.00 

13.50 
13.00 
12 . 50 [ 
12.00 

11.50 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 


H _I 

6085.45 
6118.12 
6151.04 
6184.23 
6217,70 

6251.45 

6285.47 

6319.77 

6354.35 ' 
6389.21 

6424.35 

6459.77 

6495.47 

6531.46 
6567.73 
6604.29 

6641.13 ' 
6678.26 ' 
6715.67 
6753.38 

6791.37 

6829.65 
6868 . 22 ] 
6907.09 
6946 . 24 ] 
6985.69 
7025.44 
7065.48 
7105.82 
7146 . 45 ] 

7187.38 
7228.61 

7270.14 
7311 . 97 ' 
7354.10 
7396 . 53 ] 

7439 . 27 ] 
7482 . 31 ] 

7525.66 
7569.31 

7613 . 27 ] 
7657.57 
7702.23 
7747.25 
7792 . 63 ] 
7838.37 
7884.46 
7930 . 92 ! 
7977.73 
8024.89 
8072.41 
8120.28 
8168.51 


28.02 
27.78 
27.53 
27.29 
27.05 
26.811 

26.56 

26.32 
26.07 ~ 

25.82 

25.57 

25.32 
25.07 " 

24.82 

24.57 
24.31 
24.06 
23.80 
23.55 
23.29 
23.03 
22.77 
22.52 
22.25 
21.99 
21.73 
21.47 
21.21 
20.94 
20.68 
20.41 
20.15 
19.88 
19.61 
19.35 
19.08 
18.81 
18.54 
18.27 
18.00 
17.73 
17.46 
17.19 

16.93 
16.67 

16.41 
16.16 
15.91 
15 . 6 & 

15.41 
15.17 

14.93 
14.69 


1739911.29 
1739920.87 
1739929.39 

1739936.97 

1739943.72 

1739949.73 
1739955.09 
1739959.86 
1739964 . ll | 
1739967.90 
1739971 . 27 ] 

1739974.29 

1739976.97 
1739979.37 
1739981.51 
1739983.42 
1739985 .lTT 

1739986.65 
1739988.01 
1739989.23 
1739990 . 32 ] 

1739991.29 
1739992 .nT 
1739992.95 

1739993.65 
1739994 . 27 |~ 
1739994 . 83 ] 
1739995 . 34 ] 
1739995.79 
1739996 . 20 ] 
1739996 . 56 ] 
1739996 . 89 ] 
1739997 . 18 ] 
1739997,45 

1739997.68 

1739997.90 
1739998.09 
1739998.27 
1739998 . 42 [ 
1739998.56 

1739998.69 

1739998 .TIT 

1739998.91 
1739999.01 
1739999709 ^ 
1739999.17 

1739999724] 

1739999.30 

1739999.36 

1739999.41 

1739999.46 

1739999.50 

1739999.54 


K __ 

0.001527 

0.001362 

0,001215 

0.001085 

0.000969 

0.000865 

0.000773 

0.000691 

0.000618 

0,000552 

0.000494 

0.000443 

0.000396 

0 . 000355 ] 

0 , 000318 ] 

0 . 000285 ] 

0 . 000256 ] 

0 . 000230 ] 

0 . 000206 ] 

0 . 000185 ] 

0 . 000167 ] 
0 . 000150 ] 
0 . 000135 ] 
0.00012l' 
0 . 000109 ] 
0 . 000099 ] 
0 . 000089 ] 
0 . 000080 ] 
0 . 000072 ] 
0 . 000065 ] 
0 . 000059 ] 
0 . 000054 ] 
0 . 000049 ] 
0 . 000044 ] 

0,000040] 

0 . 000036 ] 
0 . 000033 ] 
0 . 000030 ] 
0 . 000027 ] 
0 . 000025 ] 
0 . 000023 ] 

0 000021] 
0! 000019] 
0^000017] 
0 . 000016 ] 
0 . 000014 ^ 

0.000013] 
0,000012 
0.000011 
0,000010 
0.000009 
0.000009 
0.000008 


L 

0.000000 
0.000000 
0 . 000000 ] 
0,000000 
0,000000 
0.000000 
0,000000 
0 . 000000 ] 
0,000000 
0.000000 
0.0000001 
o.oooooo] 
0.000000 
0.000000 
0 .oooooo] 

0.000000] 

0.000000 

0.000000] 

0.oooooo] 
0.000000] 
0.000000] 
0.000000] 
0.oooooo] 
0.000000] 
o.oooooo] 
0 .oooooo] 
0.000000] 
o.oooooo] 
0.000000] 
o.oooooo] 
o.oooooo] 
o.oooooo] 
0.000000] 
0 .ooooo"o] 
o.oooooo] 
0.000000] 
0 .oopooT] 
o.ooooo~o] 
0 !oooooo] 
o.oooooo] 
0 oooooo] 
o!oooooo] 
o!ooooom 
o.oooooo^ 
0,oooooo 
o.oooooo 
0,000000 
0.000000 
o.oooooo 
0.000000 
0.000000 
0.000000 
0.000000 
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l / AE 448 data 



g i 

H 1 

— i—r 

J 

K 

L 

213 

11.00 

8217.09 

14.45 

1739999.58 

0.000007 

0.000000 

214 

11.00 

8266.02 

14.22 

1739999.61 

0.000007 

0.000000 

Jh a A A A A 

215 

11.00 

8315.30 

13.99 

1739999.64 

0.000006 

0 . 000000 

2 16 

11.00 

8364.93 

13.77 

1739999.67 

0.000006 

0.000000 

2 17 

11.00 

8414.92 

13.54 

1739999.70 

0.000005 

0 . 000000 

6 A r 

2 18 

11.00 

8465.25 

13.32 

1739999.72 

0.000005 0.000000 

2 18 

11.00 

8515.94 

13.10 

1739999.74 

0.000004 0.000000 

2 2 0 

11.00 

8566.97 

12.88 

1739999.76 

0.000004 

0.000000 

A» At V 

2 21 

11.00 

8618.35 

12.67 

1739999.78 

0.000004 

0.000000 

At As 

2 2 2 

11.00 

8670.08 

12.46 

1739999.79 

0.000004 

0.000000 

M A* M 

223 

11.00 

8722.16 

12.25 

1739999.81 

0.000003 

0.000000 

2 2 4 

11.00 

8774.59 

12.04 

1739999.82 

0.000003 

0.000000 

■A A» A 

2 2 5 

11.00 

8827.37 

11.84 

1739999.83 

0.000003 

0.000000 

A At tt 

2 2 6 

11.00 

8880.50 

11.64 

1739999.85 

0.000003 

0.000000 

2 2 7 

11.00 

8933.98 

11.44 

1739999.86 

0.000002 

0.000000 

2 2 8 

11.00 

8987.81 

11.24 

1739999.87 

0.000002 

0.000000 

A* At W 

2 2 8 

11.00 

9041.98 

11.05 

1739999.88 

0.000002 

0.000000 

2 30 

11.00 

9096.51 

10.86 

1739999.88 

0.000002 

0.000000 

2 31 

11.00 

9151.38 

10.67 

1739999.89 

0.000002 

0.000000 

2 32 

11.00 

9206.61 

10.48 

1739999.90 

0.000002 

0.000000 

« ^ ™ 

2 33 

11.00 

9262.19 

10.29 

1739999.90 

0.000002 

0.000000 

a* aJ m0 

2 3 4 

11.00 

9318.12 

10.11 

1739999.91 

0.000002 

0.000000 

M W A 

2 35 

11.00 

9374.40 

9.93 

1739999.92 

0.000001 

0.000000 

At At V 

2 3 6 

11.00 

9431.03 

9.75 

1739999.92 

0.000001 

0.000000 

2 37 

11.00 

9488.02 

9.58 

1739999.93 

0.000001 

0.000000 

•A A# r 

2 3 8 

11.00 

9545.36 

9.40 

1739999.93 

0.000001 

0.000000 

239 

11.00 

9603.05 

9.23 

1739999.94 

0.000001 

0.000000 

2 4 0 

11.00 

9661.10 

9.06 

1739999.94 

0.000001 

0.000000 

2 41 

11.00 

9719.51 

8.89 

1739999.94 

0.000001 

0.000000 

2 42 

11.00 

9778.28 

8.73 

1739999.95 

O.OOOOOlJ 

0.000000 

2 43 

11.00 

9837.40 

8.56 

1739999.95 

0.000001 

0.000000 

2 44 

11.00 

9896.88 

8.40 

1739999.95 

0.000001 

0.000000 

2 45 

11.00 

9956.72 

8.24 

1739999.95 

0.000001 

0.000000 

2 4 6 

11.00 

10016.92 

8.09 

1739999.96 

0.000001 

0.000000 

2 47 

11.00 

10077.49 

7.93 

1739999.96 

0.000001 

0.000000 

2 48 

11.00 

10138.42 

7.78 

1739999.96 

0.000001 

0.000000 

2 4 9 

11.00 

10199.72 

7.63 

1739999.96 

0.000001 

0.000000 

2 50 

11.00 

10261.38 

7.48 

1739999.97 

0.000001 

0.000000 

251 

11.00 

10323.41 

7.33 

1739999.97 

0.000001 

0.000000 

2 52 

11.00 

10385.81 

7.18 

1739999.97 

0.000001 

0.000000 

2 53 

11.00 

10448.58 

7.04 

1739999.97 

0.000001 

0.000000 

2 5 4 

11.00 

10511.73 

6.90 

1739999.97 

0.000000 

0.000000 

255 

11.00 

10575.24 

6.76 

1739999.97 

0.000000 

0.000000 

2 5 6 

11.00 

10639.14 

6.62 

1739999.97 

0.000000 

0.000000 

2 57 

11.00 

10703.41 

6.48 

1739999.98 

0.000000 

0.000000 

258 

11.00 

10768.06 

6.35 

1739999.98 

0.000000 

0.000000 

2 5 9 

11.00 

10833.10 

6.22 

1739999.98 

0.000000 

0.000000 

2 60 

11 . oc 

10898.52 

6.09 

1739999.98 

0.000000 

0.000000 

2 61 

11 .oc 

10964.32 

5.96 

1739999.98 

0.000000 

0.000000 

2 62 

11 . oc 

11030.52 

5.82 

1739999.98 

0.000000 

0.000000 

2 63 

11.oc 

11097 . 1 C 

5 . 7 C 

1739999.96 

0.000000 

0.000000 

.2 6 4 

11 .oc 

11164 . 0 " 

5.56 

1739999.96 

0 .oooooc 

0.000000 

teii 

11.oc 

11231.44 

5.46 

1739999.96 

0 . 00 D 00 C 

0.000000 
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l / AE 448 data 


—r 

H 1 

— 1 —[ 

J 

K 

L 

11,00 

11299.21 

5.33 

1739999.98 

0.000000 

0.000000 

11.00 

11367.38 

5.22 

1739999.98 

0.000000 

0.000000 

11.00 

11435.95 

5.10 

1739999.99 

0.000000 

0.000000 

11.00 

11504.92 

4.98 

1739999.99 

0.000000 

0.000000 

11.00 

11574.30 

4.87 

1739999.99 

0.000000 

0.000000 

11.00 

11644.10 

4.75 

1739999.99 

0.000000 

0.000000 

11.00 

11714.30 

4.64 

1739999.99 

0.000000 

0.000000 

11.00 

11784.92 
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Figure 3.0.2: Plot of X-Position vs. Time 
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Figure 3.0.3: Plot of Y-Velocity vs. Time 
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Figure 3.0.4: Plot of Y-Position vs. Time 
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Figure 3.0.5: Plot of Weight vs. Time 
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Figure 3.0.6: Plot of Inclination vs. Time 
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Figure 3.0.9: Plot of Thrust vs. Time 
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Figure 3.0.10: Plot of Thrust/Weight Ratio vs. Time 
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Figure 3.0.11: Plot of Y-Position vs. X-Position 
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4.0 Orbital Analysis 


The NASA ground rules for the Shuttle II specified two orbits into which the orbiter must be 
capable of launching: a 270 nautical mile orbit due east of Kennedy Space Center, and a 150 
nautical mile polar orbit To verify the feasibility of the proposed ASTS configuration, it 
was necessary to determine the types and magnitudes of the orbital maneuvers that would 
be necessary for orbit inserton. As a starting point, the output data for the final position and 
velocity from the launch simulation program was used. This data places the orbiter at an 
initial altitude of 78.1 nautical miles with a velocity of 25,962 feet per second. 


4.1 The Insertion Orbit 

The initial position and velocity of the orbiter place the orbiter originally at the perigee of 
an elliptical orbit. The semi major axis of this orbit, ai, was calculated using the Vis-Viva 

equation: . ^ • 

* <: 

where ri is the sum of the altitude of the orbit and the radius of the earth (20925672.57 feet). 

The eccentricity of this first orbit, and the altitude and velocity at the apogee of this orbit 
were calculated using simple geometric relationships. The pertinent data on the first orbit 

is as follows: 



ei = .02473 

ai = 21,943,244 feet 

Hap = 22,485371 feet 

Viap = 24,709 feet/second 

4.2 The 270 Nautical Mile Orbit ~ 

In order to provide a minimum energy transfer between the insertion orbit and the 270 nmi 
orbit, an elliptical transfer orbit was used (see figure 4.2). 

The 270 nmi orbit has a radius of a2 = 22,567,272 ft. The semi major axis of the transfer 
ellipse, at, has a value of (ai +a2)/2 = 22,255,258 ft. 

The velocities at periapsis and apoapsis of the transfer ellipse, as well as the velocity for 
the circular orbit were calculated using the Vis-Viva equation again. 








The following values were calculated 
insertion into a 270 nmi orbit: 


for the velocities and the velocity changes needed for 


Viap 

Vti 

dVi 


24,709 ft/sec 
24,890 ft/sec 

181 ft/sec 


Vt2 

V270 


24,799 ft/sec 
24,975 ft/sec 


dV 2 


176 ft/sec 


% 

The total impulsice change in velocity for this orbital maneuver is 351 ft/sec. 



4.3 The 150 Nautical Mile Orbit 

The path of the elliptical orbit transfer orbit used for this maneuver is shown in figure 4.3. 
The velocity and radius of the 150 nmi orbit are found to be 


ri50 

Vi50 


21837672ft 
25389 ft/sec 


Following the sane procedure as before in order to obtain the velocity changes necessary 
yields the following information: 


Viap 

Vti 

dVi 

Vt2 

Vi50 


24709 ft/sec 
24678 ft/sec 

-31 ft/sec 

25419 ft/sec 
25389 ft/sec 


dV 2 


-30 ft/sec 


The total impulsive change in velocity for this orbital transfer is -61 ft/sec. These two 
velocity changes would have to be made in a direction opposite to the direction of travel of 

the shuttle. 
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Figure 4.2.1:270 Nautical Mite Orbital Transfer 



Figure 4.3.1:150 Nautical Mile Polar Orbital Transfer 
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5.0 Stability and Control Analysis 

Stability and Control Analysis for the Advanced Space Transportation System was 
performed using geometric Shuttle data and empirical and graphical methods given by 
Etkin, Roskam and Smetana. The geometric parameters used to perform the analysis 
included physical data for the wing, the body, and the vertical tail. 

In order to determine stability and control derivatives, the geometric data for the Shuttle II 
must be used to determine aerodynamic parameters such as aspect ratio, wing taper ratio, 
and the mean aerodynamic chord. 

’A' is the aspect ratio and is defined to be equal to b^S; b is the wing span and S is the wing 
planform area. 'S' can be easily determined using the following formula: S=(c r +ct)(b/2), 
where c r is the wing root chord and ct is the wing tip chord. The ratio of the wing tip chord to 
the wing root chord is in itself and important physical parameter known simply as the 
taper ratio, X. The value of the mean aerodynamic chord (the wing chord which passes 
through the centroid of either half-wing) can be found by integrating the equation for the 
square of the wing chord: c 2 = {c r - [(c r -ct)/(b/2)]y) 2 over the span of half a wing and 
multiplying by 2/S. All of these aerodynamic parameters have been calculated and are 
foundm^tiie^ExecutiveSummaryJn^Tables^2^and£^L_ 

5.1 Determination of Clu 

CLa is an extremely important stability derivative describing the change in lift of an 
aircraft due to variation in angle of attack. 

Roskam gives the following formula for determining Clo: 

CLa = Kwh^Low 

where K w b =1 - .25(d/b) 2 + 0.25(d/b) 

and Cl«w = 2jtA/{2+[(A 2 B 2 /K 2 )(l+tan 2 Ac/2/B 2 )+4] 1/2 ) 

In these equations: d = body diameter at root chord = 30 ft. 

b = wing span = 127.7 ft. 

A = wing aspect ratio = 2.8 

P = compressibility parameter = (1 - M 2 ) 172 (for subsonic flight) 

K = ratio of actual wing section left curve slope to 2 jc 
Ac/ 2 = sweep of the wing half-chord line = 34° 

Since these equations are useful only for Mach numbers less than the critical Mach 
number, K is assumed to be approximately 1. Hence, the following values for Cl<x are 

obtained: 


M range jj2 

0.0 - 0.3 1.0 

0.3 - 0.6 0.775 

0 .6-Mcr 0.415 


(/rad) 

2.97 

3.09 

3.33 


It is apparent that these values alone for the lift-curve slope (although correct) cannot be 
used to represent Clu for the ASTS. For this reason, methods of estimating Cl« at higher 
Mach numbers and varying Reynolds numbers (from Dynamics of Flight, by Etkin) were 

employed. 





Etkin presents a graphical method by which the lift-curve slope of an aircraft may be 
obtained as a function of its geometrical parameters. Before incorporating Etkin s method, 
values for t/c and r (defined in Figure 5.1.1) were assumed to be 0.06 and 5 , respectively. 

Using Figure 5.1.1 the following section lift curve slopes were obtained: 


Ba 

lailflZlaiioT 

iailoi/radl 

10^ 

0.87 

5.74 

10 7 

0.93 

6.14 

10 8 

0.95 

6.27 


It follows that these values for the actual lift curve slope are good approximations for 
transition both at the leading edge and at c/2. 

These values of (ai)o were then used in conjunction with Figure 5.1.2 to estimate the 
overall aircraft lift curve slope. Also, Etkin defines K as [(aiW2n]p and Ap as the 

sworn narameter eaual to tan' 1 (tanAcM/P). 


Table 5.1.1. CLa of ASTS for M < Mcr 


M Range 


0.0 - 0.3 


0.0 - 0.3 


0.0-0.3 


.3 - 0.6 


0.3 - 0.6 




0.6 - Mcr 


0.6 - Mcr 


0.6 - Mcr 


.88 


0.88 


.88 


0.644 


0.644 


0.644 


Ab (°) 


42 


42 


42 


46 


46 


46 


54 


54 


54 


Rn 



K 


0.914 


0.977 


0.998 


0.804 


0.860 


0.878 


0.589 


.629 


0.643 


PA/k 

PClq/K </°) 

CLa (/rad) 

3.06 

0.051 

2.67 

2.81 

0.050 

2.80 

2.805 

0.049 

2.80 

3.06 

0.049 

2.57 

2.87 

0.048 

2.69 

2.81 

0.048 

2.74 1 

3.06 

0.0451 

2.36 

2.871 

0.043 

2.41 

2.80 

0.043 

2.46 


Note that these values of CLa are slightly higher than those estimated using the method 
outlined by Roskam. 

Etkin presents a system for predicting Cl« for Mach numbers greater than the critical 
Mach number using the graphical method shown in Figure 5.1.3 which is a function o , 
M, X, and Ac/2. Results using this figure are listed in Table 5.1.2. 
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Table 5.1.2. CLa of ASTS for M > Mcr 



Prediction ofCLa for Mach numbers peater than approximately 2.8 was not possible due to 
the lack of a proper estimation technique. 



5.2 Determination of the Neutral Point, hn 


The overall aircraft neutral point, hn, is defined as the point wheretnepucmng.H 
invariant with angle of attack. It is further defined geometrically by Figure 5.2.1. 
Roskam gives the following relationship for determining the neutral point . 

Tin = hnw + Ahnn (for tail*less configurations) 


where hnw = ki[(Xac’/cr)*k2] 

and AhnB - - (dM/da)/(qScCLaw) 

where dM/da = (q/36.5) S Wf^QCilde/daAXi 


Note that Koskam defines both the neutral point, hn, and the aerodynamic center, Xac/c, as 
the same quantity. 


Using the relationships for AhnB, it was found that the body had little or no effect on the 
overall aircraft neutral point Therefore, determination of the neutral point is reduced 

the following form: 


hn = ki[(Xac'/c r )-k 2 ] 

Xac'/cr ki and ko are found using Figures 5.2.2, 5.2.3, and 5.2.4, respectively. Using 
these figures together with known geometric data and the relationship for hn, the following 

results were obtained. 



Table 53.1. Neutral Point of the ASTS for M < Mcr 
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These values for hn show that as the Mach is increased, the neutral point begins to move 
gradually aft. This finding agrees with most aerodynamic data on similar aircr 


It is important to determine the characteristics of a flight vehicle in as many of lts 
regimes as possible. Therefore, methods given by Etkin were used to find values of h 

Mach numbers greater than the critical Mach number. 


Etkin gives a graphical method shown in Figure 5.2.5 for determining hn as a function of 
A, M ,X, and Ad2- Using these variables and Figure 5.2.5, the following data was obtained. 



Table 5A1 Neutral Point of the ASTS for 0 < M < 2Ji 



Again, these numbers show that (as expected) the neutral point starts out relatively close to 
the quarter-chord for low subsonic Mach numbers and moves aft approaching the ha 
chord at higher supersonic Mach numbers. 


5.3 Determination of Cmu 

C mo is an extremely important stability derivative which dictates whether or not a 
configuration has positive or negative static stability- Cmo is determined through the 

following equation: 


Cmo = CLa (h * hn) 


where h represents the non-dimensional location of the aircraft 

configuration to be statically stable, Cmo must be negative; this requires that the quantity 
(h - hn) must be negative. Xcg was determined by considering the location and weight of 
the components of the ASTS. The distance between the center of gravity and the leading 
edge of the mean aerodynamic chord (non-dimensionalized by dividing by the mean 
aerodynamic chord) is defined as h. Upon completion of these calculations, it is found that 

h = 0.136. 


Hence, _Cmcc = Clo (0.136 * hn) 

5.4 Clo of Vertical Tail 
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The lift curve slope of the vertical tail is determined by the geometric characteristics of the 
vertical tail and is influenced by the change in flow characteristics due to the presence of 
the wing and body combination. Figure 5.4.1 shows the relationship between (ap)B/ai for a 
vertical tail on a body of circular cross section. Using a value of 1 for D/h, (aF)B/ai is 
approximately 2.9. 


Values otfcSfor the vertical tail are determined using the same procedure as before for the 
wing. 'When-this value is obtained it is multiplied by the factor (ap)B/ai to correct for its 
location ufthe flow field relative the the wing and body. It is assumed that the vertical tail 
has approximately the same values of t/c and x as the wing so that both the wing and 
vertical tail have the same values for (ai)o. Using the geometric parameters for the 
vertical tail and Figures 5.1.2 and 5.1.3 the following is found. 



Table 54.1. Corrected Clo values for the ASTS Vertical Tail 



Determination of Clo corrected for the vertical tail beyond a Mach number of 
approximately 1.6 was not possible. 


5.5 Determination of C| P 

Clo is the change in rolling moment due to a variation in rolling velocity. Ci p is often 
referred to as the rolling moment damping derivative. Etkin provides a graphical method 
by which Ci p can be calculated based on the following parameters: P, K r A,^»a nd Ae (which 
is defined as tan* 1 [(l/pltanAc/J). Etkin’s method is presented here in Figure 5.5.1. Table 
5.5.1 contains the values of Ci p found using Etkin s method. 
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Table 5.5.1. Ci„ values ofASTS for M < Mcr 


M range 


0.0 - 0.3 


0.0 - 0.3 


0.0 - 0.3 


0.3 - 0.6 


0.3 - 0.6 


0.6 - Mcr 


0.6 - Mcr 


0.6 - Mcr 


1.0 


1.0 


.88 


0.88 


0.644 


0.644 


0.644 


Ae(°) 


42 


42 


42 


46 


46 


46 


54 


54 


54 


Rn 


K 

Cio (/rad) 

0.914 

-0.219 

0.977 

-0.225 

0.998 

-0.230 

0.804 

-0.215 

0.860 

-0.225 

0.878 

-0.229 

0.5891 

-0.201 

0.629 

-0.210 

0.643 

-0.215 


Cip does not vary greatly with Mach number, and is always negative. 

5.6 Determination of Cip 

CiB represents the change in rolling moment due to a change in sideslip angle, p. Cip is 
primarily governed by the wing dihedral angle, the body and the vertical tail. For this 
reason, Cip is often referred to as the effective dihedral derivative. Eltkin presents a 
graphical method of determining Cip based on the effects of the wing, the body and the 
vertical tail. This method utilizes Figure 5.6.1. Using this figure, it was found that 

Cip w = *.32 Cl. 

The effect of the body on C]p is found using the equation: 

(AQp)B * 1-2 (A) 1/2 Zw (h+wl/b 2 

where Zw = vertical distance of c/4 point below fuselage centerline = 11.0 ft. 
h s average fuselage height at wing root * 28.0 ft. 
w = average fuselage width at wing root = 30.0 ft. 

Using this information: (ACip)B = 0.08. 

The effect of the vertical tail on C]p is found using the equation. 


(ACip)p * -ap (1 * dc/dp)(Vp/V)2 Sp Zp/(S b) 

where Zp = vertical distance between vertical tail AC and aircraft CG = 41.0 ft. 
Sp = Vertical Tail planform area = 603 ft 2 . 

Using this equation under the assumption that da/dp is negligible and Vf/V is 
approximately equal to 1: (ACip)F = -0.033 ap 

Combining components: Cip* -0.32 Cl + 0.08 - 0.033ap 

This equation holds for .ill Mach numbers below the critical Mach number. 





5.7 Determination of Ci r 

Clr represents the change in rolling moment due to yawing velocity, or simply rolling 
moment due to yaw. Etkin divides Ci r into two components: Ci r due to the wing and Q r due 

to the vertical tail. Using Figure 5.7.1, (Ci r ) B = 0-23 C L- 

The contribution of the vertical tail on Ci r is shown through the following formula: 


(Qr )F = S F Z F fc +(Vb)]/(S b) 

where If * distance between vertical tail AC and the aircraft CG = 48.0 ft. 
Hence, (Q r )F ■ 0.079 aF 

Combining components: Ci r = 0.23 Cl + 0.079 aF 

This equation holds for all Mach numbers below the critical Mach number 



5.8 Determination of C n p 


C n B is defined to represent the change in yawing moment due to sideslip angle, (i. C n p is a 
very important stability derivative since it alone determines whether or not an aircraft has 

_ _ ... _*1 *j r* * ho nncirivp 


positive weathercock stability; for positive weathercock stability, C n p must be positive. 

Figure 5.8.1 shows the variation of C n p/C L 2 as a fimction of wing aspect ratio, A, and wing 
quarter-chord sweep angle, A. 

Using Figure 5.8.1, C n p = 0.07 Cl 2 


5.9 Determination of C n p 

C np is a cross derivative defined as the change in yawing moment due to rolling velocity. 
Etkin presents the following equation for determining Cn p : 

Cnn = KACnoh/CiJ Cl + KAC np )2/(CD0)a] (CD0)a 


where (ACnp)l/CL and (ACnp)2/(CD0)® are found f¥om 


5.9.1 


From Figure 5.9.1, (ACnp)]/CL = -0.17 and (ACnpl^CDok* 

3 S 

Cnp = -0.17 Cl + 19 (Cdo)cl 


= 19 , so that C np may be written 


5.10 Determination of Cnr 

Cnr is defined to be the change in yawing moment due to yawing velocity and is often 
referred to as the damping in yaw derivative. Etkin presents the following equation for 

determining C n r: 

Cnr = [(AC nr )i/CL 2 ]C L 2 + [(ACnr)2/ClX)]CD0 
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where [(ACnrJi/CL 2 ] and KACnr)2/CD0] are found from Figure 5.10.1 

Figure 5.10.1 reveals that [(ACnrh/CL 2 ] = -0.02 and [(ACnr^CDol = -0.48, so that C nr may be 

m 

written as 

Cm- = -0-02 Cl 2 - 0.48 CdO 

5.11 Determination of CMq 

C mq is defined as the variation in pitching moment due to pitch velocity or simply as the 
damping in pitch. Smetana gives the following relationship for Cmq : 

Cmq =-2x' lx’I Clo/c 2 

where x' = distance from CG to wing quarter-chord (positive for forward CG) = -6.97 ft. 

After substituting the proper geometric and stability parameters, C m q can be written as 
Cmq = 0.0372 Cjxt 

5.12 Determination of Ciq 

CLq is defined as the change in lift coefficient due to pitch velocity. Smetana gives the 
following relationship for CLq : 

CLq = 2 x Cl<x/ c 

Upon substitution of known geometric parameters, CLq is written as 
CLq = -0.2727 Clo 

5.13 Determination of Cd<x 

Cdo is defined as the change in drag coefficient due to a change in the angle of attack. 
Smetana gives the following relationship for CDa: 

CDa = (dCDo/da) + [2 Cl/(*A)] Clo 

Since (dCD 0 /da) is usually small, Cdo may be written simply as [2 Cl/0eA)] c La- 

Therefore, CDa may be expressed as 
Cdo = 0.2274 Cl Clo- 

5.14 Determination of Cyp 

C y B is defined as the change in side force,y, due to a change in sideslip angle, p. Smetana 
breaks Cyp up into 3 components: that of the wing, body, and vertical tail. 

(Cyp)w = Cl 2 (6 tan A sin A)/[A(A + 4 cosA)] 

Upon substitution, (Cyp) w = 0.2232 Cl 2 
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(C y p)B = -ki (Clo)b (Sb/S) 

where Sb = body reference area = 3733 ft* 

ki = interference factor determined from Figure 5.14.1 


Figure 5.14.1 gives ki as a function of Zw/(d/2) where Zw is the distance froin 0 y 
centerline to quarter-chord point of exposed wing root chord (positive for quarter-chord 
point below body centerline) and d is the max body height at wmg/body intersection. 


Zw = 11.0 ft 
d = 25.0 ft 


From Figure 5.14.1, ki = 1.44 

Smetana gives the following values of (Cl<x)b : 

(Cl<i) circular body = 0.0525 rad' 1 
(Clo) rectangular body = 0.1253 rad' 1 

Since the ASTS configuration is a combination of these two shapes, let (Clo)b = °- 08 rad 1 

Now that all values needed are known, (Cyp)fi = -0.074 

% 

(Cyp)F = ' k ( C Lo)F d + do/d P> S F /S 

k is found from Figure 5.14.2 and is a function of IW(2ri)]. 

by = vertical tail span = 29.0 ft. 

ri = fu s e lag e diameter in tail area = 30 ft. 

From Figure 5.14.2, k = 0.75 

If it is assumed that do/dB is negligible then (Cyp)p may be expressed as 

(Cyp)F = '°- 078 a F 

Combining all components: Cyp * 0.2232 Cl^ • 0.074 - 0.078 ap 

5.15 Determination of C yp 

Cyp represents the change in side force due to rolling velocity. Smetana expresses C yp as 
the sum of the contributions due to the wing and the vertical tail. 

(Cyp ) w = Cl (t(A + cosA) tanMA + 4 cosA)] + 1/A) 

Substitution into (Cyp )w yields (Cyp )w = 0.91 Cl 


(Cyp )f = -2 ap Sp Z F /(S b) 

All of these parameters have been defined previously so that (Cyp )F = -0.067 aF 
Combining these two: Cyp = 0.91 Cl - 0.067 aF 
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5.16 Determination of Cyr 

Cyr is defined as the change in side force due to yawing velocity. Smetana expresses Cyr 
asthe sum of the contributions due to the wing and the vertical tail. 


(Cyr)w = 0.143 Cl • 0.05 

(Cyr)P = -2 If (Cyr)F/b = -0.752 (Cyr)p 


So that, Cyr = 0.143 Cl- 0.05 - 0.752 (Cyr>F 

5.17 Dynamic Stability Matrices 


Once the proper stability derivatives have been determined, it is possible to use them to 
investigate the dynamic stability of the configuration. Presented here are the longitudinal 
and lateral stability matrices corresponding to the non-dimensional, controls fixed, 
linearized equations of motion for a rigid aircraft. 


Longitudin al Equations: 

Dx = Ax where x = (u a q 0) T 

A(l,l) = (Cxu + 2 Clo tan©o)/2n 
A(l,2) = CxoAi 
A(l,3) = 0.0 
A(l,4) = -C\jyt2\i 

A(2,l) = (CZu * 2 CL0V(2p. - Czo) 
A(2,2) = Cz<V(2|X - Czo) 

A(2,3) = (2p. + CZqV(2p. - Czo) 

A(2,4) = -Clo tan©o/(2p - Cz«) 

A(3,l) = (C M u + Cmo A(2,l))/iB 
A(3,2) = (CMa + C m a A(2,2)ViB 
A(3,3) = (CMq + Cma A(2,3))/iB 
A(3,4) — Cma A(2,4)/iB 
A(4,l) = 0.0 
A(4,2) = 0.0 
A(4^) = 1.0 
A(4,4) = 0.0 

where p = m/(p S c/2) 

iB - iyy/tp S (c/2) 3 ] 

CXa = CLO * CDa 
Cza = -Clo * Cdo 
C zu = - (M 2 /l-M 2 ) C lo 
CZ q = -CLq 
CDm = CDu/M 

Cxu = -2 (Cdo + Clo tan0o) - M Com 

Czo = 0.0 



62 



Dy = By where y = (B p r <bV 

B(l,l) = Cyp/2n 
B(l,2) = Cyp/2M- 

B(l,3) = -(2p - Cyr)/2|i 
B(l,4) = Clo^P 
B(2,l) = (ic Cip + iE C n pVA 
B(2,2) = ttc Ci p + i£ C np )/A 
B(2,3) = (ic Or + iE C nr VA 
B(2,4) = 0.0 

B(3,1) = Ga C n p + iE QpVA 
B(3,2) = Ga C np + iE QpVA 
B(3,3) = Ga Qir + iE QrVA 
B(3,4) = 0.0 

B(4,l) = 0.0 
B(4,2) = 1.0 
B(4,3) = tanQo 
B(4,4) = 0.0 

where p. = m/(p S h/2) 

A = iAiC-iE 2 
iA = Ixx/[p S (b/2) 3 ] 
ic = Izz/[p S (h/2) 3 ] 
ig = Ixz/[p S (h/2) 3 ] 

Once the matrices A and B are found, their Eigen Values can be determined. These Eigen 
values reflect whether or not the system is stable. For a system to asymptotically stable, all 

of its Eigen values must have negative real parts. 

The matrices A and B were calculated using a computer program written in FORTRAN. 
The program and results of the investigation of the dynamic stability for landing 
(assuming 0 O = 10° and U 0 = 300 ft/sec and L = W) are listed on pages 66 and 71. 

5.18 Determination of C|§ a 

QSo describes the variation in rolling moment with change in aileron deflection; often, 
this stability derivative is referred to as the aileron power. Smetana gives the following 

equation for Ci5a- 

Ab 

Q6a = [2CLax/(Sb)]\ cydy 

where c is the equation for the wing chord as a function of wing span, y. 


c = c r - [(c r -c t V(b/2)] y = 73.125 - 0.863 y 

and a and b represent the spanwise distance for the body centerline to the inner and outer 
aileron points, respectively. 







a = 27.0ft b = 58.5 ft 

x is a correction factor which is a function of the aileron chord to wing chord ratio, x is 
determined from Figure 5.18.1. Since the ratio of aileron chord to wing chord is not 
constant over the span of the aileron, the average value was used to determine x. 

(c a /c w )avg = 0.33 fcs0.57 

Upon integration and substitution of the proper values, 

Ciso = 0.071 Clo. 

5.19 Determination of C|gr 

Ci5r is the variation in rolling moment coefficient with change in rudder deflection. 
Smetana gives the following formula for estimating C]5r : 


CiSr = ap x Sp Zp/(S b) 

Here, x is a function of the rudderarea to vertical tail area ratio 
Sr = 157.5/fpj) Sp a 603 ft2 

Hence, Sj^Sp = 0.261 so that from Figure 5.19.1, x = 0.48. 
Substitution into the equation for CiJr results in 


Cwr s 0.0158 ap. 

5.20 Determination of C n sr 

Cn5r is the variation in yawing moment coefficient with a change in rudder deflection and 
is often referred to as the rudder power. Smetana gives the following relationship for C n 6r 

C n Sr = *aF x Sf lf/(S b) 

As with Ci 5 r, the value of x is determined from Figure 5.19.1 and has a value of 0.48 so that 
C n 5r may be expressed as 

Cngr * -0.0187 ap 
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AE448 CG Locations/I 




18 


Component 





Vertical Tail 


Forward Bod’ 


Mid Bod 


Aft Bod 


Main Engines _ 


QMS __ 


RCS __ 


Control Systems 


Avionics _ 


Personnel _ 


Stores _ 


Engine Accessories 


Hydrogen Tanks_ 


Oxygen Tank _ 


TPS __ 


Totals __ 


Centroids 


Location 


116.60 


140.00 


28.90 


80.00 


128.00 


137.00 


138.00 


138.00 


48.00 


22.00 


35.00 



91.84 


c 

D 

Waiaht 

Y-Moment 

28,300 

3,299,780 

3,200 

448,000 

27,000 

780,300 

23,600 

1.888,000 

13,200 

1,689,600 

22,800 

3,123,600 

3,700 

510,600 

1,600 

220,800 

4,000 

192,000 

3,100 

68,200 

1,600 

56,000 

600 

25,200 

6, 200 

775,000 

24,000 

1,536,000 

11,000 

1,331,000 

6,100 

587,430 

180,000 

16,531,510 
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30 


DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 


PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 


QyS v W t US *CBARy UO yRHO»THETAO 
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FORMAT! ////y 1 OX ; 'NON-DIMENSIONAL STABILITY DFr’IL-<TT JEo * t 
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' LATERAL-DIRECTIONA1. I'ERIVATIV 
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' CLR 
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LL0NG=CBAR/2♦0 
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CXA=CLO-CDA 

CZA=-CLA-CDO 

S0UND=1116*4 

M=UO/SOUND 

CZU=-<(M#M)/(1♦O-MKM))#CLO 

CZQ=-CLG 

CDM--CDU/M 

CXU=- (2♦0# < CD0 + CL0#TAN(THETAO) )) 
CZAD=0.0 


-(M#CDM) 


ALONG ( 
ALONG< 
ALONG< 
ALONG( 
ALONG( 
ALONG( 
ALONG< 
ALONG( 
ALONG( 
ALONG( 

ALONG( 
ALONG( 
ALONG( 
ALONG( 


1 ? 1 ) = <CXU+2.0*CL0*TAN<THETA0>>/<2.0*MUL0NG> 

1 ?2)=CXA/(2.0#MUL0NG) 

1 ? 3)=0♦0 

1 „ 4 )=-CLO/(2♦0*MUL0NG) 

2 ?1) = (CZU-2.0*CL0)/(2♦0*MUL0NG-CZAD) 
2?2)=CZA/(2.0*MUL0NG-CZAD) 

2?3)=(2.0#MUL0NG+CZQ)/< 2 ♦0*MUL0NG-CZAD) 

2 f 4)=(-CLO*TAN(THETAO) ) / < 2.0*MUL0NG-CZAD) 

3 i) = ( CMUTCMAD*ALONG ( 2? 1) )/IBLONG 
3,2) = (CMAtCMAD#ALONG(2 ? 2))/IBLONG 
3 ? 3) = (CMQFCMAD#ALONG(2 ? 3))/IBLONG 

CMAD#ALONG(2 


4 ? 


4 ) = 


,4)/IBLONG 


0.0 

0.0 


4 ? 1 ) 

4 ? 2) 

ALONG(4 ? 3)-1 ♦ 0 

? 4 ) r - 0 » 0 


ALONG <4 - 1 


*;!<* CALCULATIONS FOR LATERAL MOTION 
LLATaWS/2.0 

MJLAT=(W/G)/(RKO#S#LLAT) 

IALAT=A/< RH0*S#LLAT#*3) 

ICLAT=C/(RHC#S*LLAT#*3) 
IELAT-E/<RH0#S#LLAT##3) 

DEL=(IALAT&ICLAT)~IELAT##2 


C 


\ 

\ 


cn 


\ 



ALAT ( 
ALAT ( 
ALAT ( 
ALAT ( 
ALAT ( 
ALAT ( 
ALAT ( 
ALAT ( 

ALAT ( 
ALAT ( 
ALAT ( 
ALAT ( 
ALAT ( 
ALAT ( 
ALAT ( 


1 y 1 ) 

1 9 2) 

1 y 3) 
1 9 4) 
1 ) 




2 y 2) 
2 y 3) 


2 y 


4) 

3 y 1 ) 
3 y 2) 
3 y 3 ) 

3 y 4) 

4 y 1) 
4 y 2) 
4 y 3) 


SLTB7 T'474) 


~CYB/ ( 2 ♦ 05KMUL.AT ) 

= CYP/ ( 2 ♦ O&MtJLAT ) 

a- ( 2 * 0>KMULAT-CYR ) / < 2. 0#MULAT ) 

«CL0/<2*0*MULAT) 

■-.( IC L A T * C L B +1E L A T * C N B) / D E L 
s (ICLAT#CLP+IELAT*CNP) /DEL 
»(ICLAT#CLR+IELAT*CNR)/DEL 
--=0 * 0 

=<IALAT#CNB+IELAT*CLB)/DEL 
=(IALAT#CNP+IELAT*CLP)/DEL 
=(IALAT#CNR+IELAT*CLR)/DEL 
= 0*0 
= 0*0 
= 1*0 

=TAN<THETA0) 

=0 


10 

c 


15 

C 

100 


101 


200 

201 

i. 


WRITE(10 y100) 

DO 10 J = 1y 4 

WRITE(10 y 200) < ALONG(J y K)y 

IF < J♦NE♦4) WRITE (10y201) 

CONTINUE 

WRITEClOy101 ) 

DO 15 J=1y 4 

WRITE(10 y 200) (ALAT(J y K)y 

IF(J♦NE * 4) WRITE (10»201) 
CONTINUE 


K=1 y 4) 


K = 1y 4) 


FORMAT <//////y 25X y ' THE STABILITY MATRIX CORRESPONDING TO',/, 

* 20X y'THE NONDIMENSIONALy CONTROLS FIXED, LINEARIZED,'y/, 

* 25Xy 'LONGITUDINAL EQUATIONS OF MOTION FOR'y/y 

* 35Xy'A RIGID AIRPLANE',//) 

FORMAT (.////9 25X y ' THE STABILITY MATRIX CORRESPONDING TO ' t /» 

* 20Xy'THE NGNDIMENSIONALy CONTROLS FIXEDy LINEARIZEDy' */ » 

* 23Xy'LATERAL EQUATIONS OF MOTION FCR'y/y 

* 35Xy'A RIGID AIRPLANE'y//) 

FORMAT<5Xy' ! 'y4<F12.3y6X)y ' 

FORMAT(5X» ' I ' y?4Xy'S') 

CLOSE(9) 

CLOSE<10) 

ST CP 
END 


: ' ) 
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GEOMETRICAL * INERTIAL & EQUILIBRIUM CHARACTERIST1CS l 


GRAVITATIONAL ACCELERATION (FT/SEC 2) 

WING AREA (FT~2) « 
WEIGHT (LB) = 
WING SPAN (FT) = 
MEAN AERODYNAMIC CHORD (FT) = 

AIRSPEED (FT/SEC) = 
AIR DENSITY (5LUGS/FT~3) * 
INITIAL THETA (RAD) * 
IYY (SLUG-FT~2) = 
IXX (SLUG-FT~2) = 
IZZ (SLUG-FT"2) = 
IXZ (SLUG-FT"2) = 

CLO (-) “ 

CDO (-) = 


32♦174000 

5818.000000 
180000.000000 
127.692000 
51.120000 


300.000000 

0.002377 

0.174500 

12490393.000000 

1686164.000000 

13399266.000000 

795369.000000 

0.418000 

0.050000 



NON-DIMENSIONAL STABILITY DERIVATIVES 4 . 


CMU 

(-) 

ss 

CMA 

(-) 

= 

CM AD 

( —) 


CMG 

(-) 

= 

CLA 

(-) 

=: 

CLO 

(-) 

= 

CI'A 

(-) 


CDU 

(-) 

=: 


0.000000 
-0.403200 
0.000000 
0.104200 
2.800000 
-0.763500 
0.266000 
0.000000 


LATERAL-DIRECTIONAL 


derivatives; 


CLB 

(/RAD) 


-0.338000 

CLP 

(/RAD) 

sr 

-0.225000 

CLR 

(/RAD) 

= 

0.775500 

CNB 

(/RAD) 


0.012000 

CNP 

(/RAD) 


-0.070 00 0 

CNR 

(/RAD) 

=: 

-0.030000 

CYB 

(/RAD) 


-0.706000 

GYP 

(/RAD) 


-0.195800 

CYR 

(/RAD) 

— 

0 ♦ J 1 M U 0 


\ 

\ 

\ 

\ 
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-0*00315899 

-0.02743687 

0.00000000 

0,00000000 


-0.05570903 

-0.74066202 

-0.04074188 

0.00000000 


THE STABILITY MATRIX CORRESPONDING TO _ 

T H F N 0 N DIM E N S10 N A L s> CONTROLS FIXE D * LINE A R1 L h. D •> 
LONGITUDINAL EQUATIONS OF MOTION FOR 

A RIGID AIRPLANE 


0.00480166 
-0,09003117 
-0.00745438 
0.00000000 


0.00000000 

1.02411888 

0.00192645 

1.00000000 


-0,0023278 


J 


0,00000000 

0.00000000 


THE STABILITY MATRIX CORRESPONDING TO 
THE NONDIMENSIONAL i CONTROLS FIXED» LINEARIZED i 

LATERAL EQUATIONS OF MOTION FOR 

A RIGID AIRPLANE 

♦ 


-0.01545018 

-0.50320992 

-0.04867205 

1.00000000 


-0.95942552 

1,69903754 

0.09279544 

0.17629303 


0.03298353 

0.00000000 

0.00000000 

0.00000000 
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LINEAR SYSTEM ANALYSIS RESULTS 
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LINEAR SYSTEK ANALYSIS RESULTS 
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t/c, thicknass/chord ratio 


Figure 5.1.1. Lift curve slopes for two-dimensional incompressible flow. 







, per degree 




A * 


Figure 5.1.2.Lift-curve slope for swept and tapered wings 

number. 


at speeds below critical Mach 
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Figure 5.1.3. Lift-curve slope for swept and tapered wings at supersonic speeds. 
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TAN A le 


Location of Lifting Surfaces. 
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o 30 w 

Figure 5.2.4. Aerodynamic Center Transformation Constant K 2 
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Figure 5 2.5. Chordwise position of the mean aerodynamic center of swept and tapered 
^ wings at high speeds expressed as a fraction of the mean aerodynamic 

chord. 
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\ t , degrees 


Figure 5.5.1. Q p for straight-tapered wings. 
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Figure 5.6.1. Cip for straight-tapered wings with no dihedral. 
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Aspect ratio, A 


Figure 5.7.1. Charts 


for estimating Ci r for subsonic incompressible flow. 
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Aspect ratio, A 








Aspect ratio, A 


5.9.1. Variation of (AC n p)l/CL and (AC np ) 2 /(CD0)a 

incompressible flow. 


with aspect ratio for subsonic 
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(AC. ) /C 






Aspect ratio, A 



Figure 5.10.1. Charts 


for estimating (AC nr )l/CL 2 and (ACnrWCDO 


f 


a £ 


A, degrees 

















.1. Values for wing-fuselage interference 





Figure 5.14.2. Values for k as a function of the ratio of vertical tail span to fuselage 

disinter in the tail region. 
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Figure 5.18.1.Values for t as a function of aileron chord to wing chord ratio. 



Figure 


5.19.1. Values for t hs a function of rudder area 


to vertical tail area ratio. 




